SherbyNet - /page/wikibook/entry.mu

<!DOCTYPE html>

<title>Large eddy simulation</title>

</head>

Large eddy simulation

</h1>

</header>

Large eddy simulation (LES) is a mathematical model for <a href="Turbulence" title="Turbulence">turbulence</a> used in <a href="Computational_fluid_dynamics" title="Computational fluid dynamics">computational fluid dynamics</a>. It was initially proposed in 1963 by <a href="Joseph_Smagorinsky" title="Joseph Smagorinsky">Joseph Smagorinsky</a> to simulate atmospheric air currents,<a href="#cite_note-Smagorinsky_1963-1">[1]</a> and first explored by Deardorff (1970).<a href="#cite_note-Deardorff_1970-2">[2]</a> LES is currently applied in a wide variety of engineering applications, including <a href="Combustion" title="Combustion">combustion</a>,<a href="#cite_note-Pitsch_2006-3">[3]</a> acoustics,<a href="#cite_note-Wagner_2007-4">[4]</a> and simulations of the atmospheric boundary layer.<a href="#cite_note-Sullivan_1994-5">[5]</a>

The simulation of turbulent flows by numerically solving the <a href="Navier%E2%80%93Stokes_equations" title="Navier–Stokes equations">Navier–Stokes equations</a> requires resolving a very wide range of time and length scales, all of which affect the flow field. Such a resolution can be achieved with <a href="Direct_numerical_simulation" title="Direct numerical simulation">direct numerical simulation</a> (DNS), but DNS is computationally expensive, and its cost prohibits simulation of practical engineering systems with complex geometry or flow configurations, such as turbulent jets, pumps, vehicles, and landing gear.

The principal idea behind LES is to reduce the computational cost by ignoring the smallest length scales, which are the most computationally expensive to resolve, via <a href="Filter_(signal_processing)" title="Filter (signal processing)">low-pass filtering</a> of the Navier–Stokes equations. Such a low-pass filtering, which can be viewed as a time- and spatial-averaging, effectively removes small-scale information from the numerical solution. This information is not irrelevant, however, and its effect on the flow field must be modelled, a task which is an active area of research for problems in which small-scales can play an important role, such as near-wall flows,<a href="#cite_note-piomelli-02-6">[6]</a><a href="#cite_note-spalart-09-7">[7]</a> reacting flows,<a href="#cite_note-Pitsch_2006-3">[3]</a> and multiphase flows.<a href="#cite_note-fox-12-8">[8]</a>

<div class="mw-heading mw-heading2"><h2 id="Filter_definition_and_properties">Filter definition and properties</h2></div>

/* start https://en.wikipedia.org/ */

.mw-parser-output .hatnote{font-style:italic}.mw-parser-output div.hatnote{padding-left:1.6em;margin-bottom:0.5em}.mw-parser-output .hatnote i{font-style:normal}.mw-parser-output .hatnote+link+.hatnote{margin-top:-0.5em}@media print{body.ns-0 .mw-parser-output .hatnote{display:none!important}}

/* end https://en.wikipedia.org/ */

</style><div role="note" class="hatnote navigation-not-searchable">Main article: <a href="Filter_(large_eddy_simulation)" title="Filter (large eddy simulation)">Filter (large eddy simulation)</a></div>

An <a href="Filter_(large_eddy_simulation)" title="Filter (large eddy simulation)">LES filter</a> can be applied to a spatial and temporal field <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \phi ({\boldsymbol {x}},t)}">

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \phi ({\boldsymbol {x}},t)}</annotation>

</semantics>

</math><img src="./f16cf44410102a86023cdc53feaa55931e106cba.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:6.6ex; height:2.843ex;" alt="{\displaystyle \phi ({\boldsymbol {x}},t)}" loading="lazy"> and perform a spatial filtering operation, a temporal filtering operation, or both. The filtered field, denoted with a bar, is defined as:<a href="#cite_note-Pope_2000-9">[9]</a><a href="#cite_note-Sagaut_2006-10">[10]</a>

<mover>

<mrow>

</mrow>

</mrow>

</mover>

</mrow>

</mrow>

</mrow>

</msubsup>

</mrow>

</mrow>

</mrow>

</msubsup>

</mrow>

</mrow>

</mrow>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\overline {\phi ({\boldsymbol {x}},t)}}=\displaystyle {\int _{-\infty }^{\infty }}\int _{-\infty }^{\infty }\phi ({\boldsymbol {r}},\tau )G({\boldsymbol {x}}-{\boldsymbol {r}},t-\tau )d\tau d{\boldsymbol {r}}}</annotation>

</semantics>

</math><img src="./4c60334c1dab18d66ea7f5c35b8de814e18ba118.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.505ex; width:46.158ex; height:6.009ex;" alt="{\displaystyle {\overline {\phi ({\boldsymbol {x}},t)}}=\displaystyle {\int _{-\infty }^{\infty }}\int _{-\infty }^{\infty }\phi ({\boldsymbol {r}},\tau )G({\boldsymbol {x}}-{\boldsymbol {r}},t-\tau )d\tau d{\boldsymbol {r}}}" loading="lazy"></dd></dl>

where <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle G}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle G}</annotation>

</semantics>

</math><img src="./f5f3c8921a3b352de45446a6789b104458c9f90b.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.827ex; height:2.176ex;" alt="{\displaystyle G}" loading="lazy"> is the filter convolution kernel. This can also be written as:

<mover>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\overline {\phi }}=G\star \phi .}</annotation>

</semantics>

</math><img src="./9a0b7f599c3da3b4507ad0fcbd74b14e8f94afc7.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:10.653ex; height:3.343ex;" alt="{\displaystyle {\overline {\phi }}=G\star \phi .}" loading="lazy"></dd></dl>

The filter kernel <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle G}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle G}</annotation>

</semantics>

</math><img src="./f5f3c8921a3b352de45446a6789b104458c9f90b.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.827ex; height:2.176ex;" alt="{\displaystyle G}" loading="lazy"> has an associated cutoff length scale <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \Delta }">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \Delta }</annotation>

</semantics>

</math><img src="./32769037c408874e1890f77554c65f39c523ebe2.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.936ex; height:2.176ex;" alt="{\displaystyle \Delta }" loading="lazy"> and cutoff time scale <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \tau _{c}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \tau _{c}}</annotation>

</semantics>

</math><img src="./ea61c9e95de7b5b55eb1ef4b3e03290002e91089.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:1.96ex; height:2.009ex;" alt="{\displaystyle \tau _{c}}" loading="lazy">. Scales smaller than these are eliminated from <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\overline {\phi }}}">

<mover>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\overline {\phi }}}</annotation>

</semantics>

</math><img src="./19839359a0b8f888d007de4891a59f29984d41e3.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:1.5ex; height:3.343ex;" alt="{\displaystyle {\overline {\phi }}}" loading="lazy">. Using the above filter definition, any field <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \phi }">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \phi }</annotation>

</semantics>

</math><img src="./72b1f30316670aee6270a28334bdf4f5072cdde4.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:1.385ex; height:2.509ex;" alt="{\displaystyle \phi }" loading="lazy"> may be split up into a filtered and sub-filtered (denoted with a prime) portion, as

<mover>

</mover>

</mrow>

<msup>

</mrow>

</msup>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \phi ={\bar {\phi }}+\phi ^{\prime }.}</annotation>

</semantics>

</math><img src="./786d4c15c126a71ddef29b0ce657c806d3d1c96b.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:11.509ex; height:2.843ex;" alt="{\displaystyle \phi ={\bar {\phi }}+\phi ^{\prime }.}" loading="lazy"></dd></dl>

It is important to note that the <a href="Filter_(large_eddy_simulation)" title="Filter (large eddy simulation)">large eddy simulation filtering operation</a> does not satisfy the properties of a <a href="Reynolds_operator" title="Reynolds operator">Reynolds operator</a>.

<div class="mw-heading mw-heading2"><h2 id="Filtered_governing_equations">Filtered governing equations</h2></div>

The governing equations of LES are obtained by filtering the <a href="Partial_differential_equations" class="mw-redirect" title="Partial differential equations">partial differential equations</a> governing the flow field <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \rho {\boldsymbol {u}}({\boldsymbol {x}},t)}">

</mrow>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \rho {\boldsymbol {u}}({\boldsymbol {x}},t)}</annotation>

</semantics>

</math><img src="./772977476318dc98fbba362bc3292848fb89bdf0.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:7.999ex; height:2.843ex;" alt="{\displaystyle \rho {\boldsymbol {u}}({\boldsymbol {x}},t)}" loading="lazy">. There are differences between the incompressible and compressible LES governing equations, which lead to the definition of a new filtering operation.

<div class="mw-heading mw-heading3"><h3 id="Incompressible_flow">Incompressible flow</h3></div>

For <a href="Incompressible_flow" title="Incompressible flow">incompressible flow</a>, the <a href="Continuity_equation#Fluid_dynamics" title="Continuity equation">continuity equation</a> and Navier–Stokes equations are filtered, yielding the filtered incompressible continuity equation,

<mfrac>

<mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\frac {\partial {\bar {u_{i}}}}{\partial x_{i}}}=0}</annotation>

</semantics>

</math><img src="./9d716c7e0ef46dfd4ec0cec3925a91b0c2920aab.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.338ex; width:8.544ex; height:5.843ex;" alt="{\displaystyle {\frac {\partial {\bar {u_{i}}}}{\partial x_{i}}}=0}" loading="lazy"></dd></dl>

and the filtered Navier–Stokes equations,

<mfrac>

<mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mrow>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mrow>

<mover>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

</mover>

</mrow>

</mrow>

<mfrac>

</mfrac>

</mrow>

<mfrac>

<mrow>

<mover>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mrow>

<mrow>

<mfrac>

<mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

</mrow>

<mfrac>

</mfrac>

</mrow>

<mfrac>

<mrow>

<mover>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\frac {\partial {\bar {u_{i}}}}{\partial t}}+{\frac {\partial }{\partial x_{j}}}\left({\overline {u_{i}u_{j}}}\right)=-{\frac {1}{\rho }}{\frac {\partial {\overline {p}}}{\partial x_{i}}}+\nu {\frac {\partial }{\partial x_{j}}}\left({\frac {\partial {\bar {u_{i}}}}{\partial x_{j}}}+{\frac {\partial {\bar {u_{j}}}}{\partial x_{i}}}\right)=-{\frac {1}{\rho }}{\frac {\partial {\overline {p}}}{\partial x_{i}}}+2\nu {\frac {\partial }{\partial x_{j}}}{\bar {S}}_{ij},}</annotation>

</semantics>

</math><img src="./6dd472db9259707a7405b2b5170b53135dd7a072.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.671ex; width:77.881ex; height:6.509ex;" alt="{\displaystyle {\frac {\partial {\bar {u_{i}}}}{\partial t}}+{\frac {\partial }{\partial x_{j}}}\left({\overline {u_{i}u_{j}}}\right)=-{\frac {1}{\rho }}{\frac {\partial {\overline {p}}}{\partial x_{i}}}+\nu {\frac {\partial }{\partial x_{j}}}\left({\frac {\partial {\bar {u_{i}}}}{\partial x_{j}}}+{\frac {\partial {\bar {u_{j}}}}{\partial x_{i}}}\right)=-{\frac {1}{\rho }}{\frac {\partial {\overline {p}}}{\partial x_{i}}}+2\nu {\frac {\partial }{\partial x_{j}}}{\bar {S}}_{ij},}" loading="lazy"></dd></dl>

where <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\bar {p}}}">

<mover>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\bar {p}}}</annotation>

</semantics>

</math><img src="./27bcf72dd4d0cc7bdee4b0ee758d4db57f2a7be6.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; margin-left: -0.089ex; width:1.449ex; height:2.343ex;" alt="{\displaystyle {\bar {p}}}" loading="lazy"> is the filtered pressure field and <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\bar {S}}_{ij}}">

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\bar {S}}_{ij}}</annotation>

</semantics>

</math><img src="./eecbfb4d27a221a88fcd9c9c2f3c04bdb5de1900.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:3.05ex; height:3.343ex;" alt="{\displaystyle {\bar {S}}_{ij}}" loading="lazy"> is the rate-of-strain tensor evaluated using the filtered velocity. The <a href="Nonlinear" class="mw-redirect" title="Nonlinear">nonlinear</a> filtered <a href="Advection" title="Advection">advection</a> term <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\overline {u_{i}u_{j}}}}">

<mover>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\overline {u_{i}u_{j}}}}</annotation>

</semantics>

</math><img src="./a850035c59b486247399b36c5ef57af0a4c88a37.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:4.484ex; height:3.009ex;" alt="{\displaystyle {\overline {u_{i}u_{j}}}}" loading="lazy"> is the chief cause of difficulty in LES modeling. It requires knowledge of the unfiltered velocity field, which is unknown, so it must be modeled. The analysis that follows illustrates the difficulty caused by the nonlinearity, namely, that it causes interaction between large and small scales, preventing separation of scales.

The filtered advection term can be split up, following Leonard (1975),<a href="#cite_note-Leonard_1974-11">[11]</a> as:

<mover>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

</mover>

</mrow>

<msub>

</mrow>

</msub>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\overline {u_{i}u_{j}}}=\tau _{ij}+{\overline {u}}_{i}{\overline {u}}_{j}}</annotation>

</semantics>

</math><img src="./bf401952ec5bf85ab8238491577419bc3473a039.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:17.515ex; height:3.009ex;" alt="{\displaystyle {\overline {u_{i}u_{j}}}=\tau _{ij}+{\overline {u}}_{i}{\overline {u}}_{j}}" loading="lazy"></dd></dl>

where <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \tau _{ij}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \tau _{ij}}</annotation>

</semantics>

</math><img src="./97d0190e7517254a655e8033badbc01c8590381c.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:2.493ex; height:2.343ex;" alt="{\displaystyle \tau _{ij}}" loading="lazy"> is the residual stress tensor, so that the filtered Navier-Stokes equations become

<mfrac>

<mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mrow>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mrow>

<mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mrow>

</mrow>

<mfrac>

</mfrac>

</mrow>

<mfrac>

<mrow>

<mover>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\frac {\partial {\bar {u_{i}}}}{\partial t}}+{\frac {\partial }{\partial x_{j}}}\left({\overline {u}}_{i}{\overline {u}}_{j}\right)=-{\frac {1}{\rho }}{\frac {\partial {\overline {p}}}{\partial x_{i}}}+2\nu {\frac {\partial }{\partial x_{j}}}{\bar {S}}_{ij}-{\frac {\partial \tau _{ij}}{\partial x_{j}}}}</annotation>

</semantics>

</math><img src="./5f388a6f625689c2a9793f3584e3fc406de88e70.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.671ex; width:49.707ex; height:6.509ex;" alt="{\displaystyle {\frac {\partial {\bar {u_{i}}}}{\partial t}}+{\frac {\partial }{\partial x_{j}}}\left({\overline {u}}_{i}{\overline {u}}_{j}\right)=-{\frac {1}{\rho }}{\frac {\partial {\overline {p}}}{\partial x_{i}}}+2\nu {\frac {\partial }{\partial x_{j}}}{\bar {S}}_{ij}-{\frac {\partial \tau _{ij}}{\partial x_{j}}}}" loading="lazy"></dd></dl>

with the residual stress tensor <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \tau _{ij}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \tau _{ij}}</annotation>

</semantics>

</math><img src="./97d0190e7517254a655e8033badbc01c8590381c.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:2.493ex; height:2.343ex;" alt="{\displaystyle \tau _{ij}}" loading="lazy"> grouping all unclosed terms. Leonard decomposed this stress tensor as <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \tau _{ij}=L_{ij}+C_{ij}+R_{ij}}">

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \tau _{ij}=L_{ij}+C_{ij}+R_{ij}}</annotation>

</semantics>

</math><img src="./e6dd698b03e8d5c4952acf4647c66c3899228ddd.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:20.713ex; height:2.843ex;" alt="{\displaystyle \tau _{ij}=L_{ij}+C_{ij}+R_{ij}}" loading="lazy"> and provided physical interpretations for each term. <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle L_{ij}={\overline {{\bar {u}}_{i}{\bar {u}}_{j}}}-{\bar {u}}_{i}{\bar {u}}_{j}}">

<msub>

</mrow>

</msub>

<mover>

<mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mrow>

</mover>

</mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle L_{ij}={\overline {{\bar {u}}_{i}{\bar {u}}_{j}}}-{\bar {u}}_{i}{\bar {u}}_{j}}</annotation>

</semantics>

</math><img src="./fc8ac43ad61212c6e48d21df43836b0708fe2ae3.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:17.851ex; height:3.343ex;" alt="{\displaystyle L_{ij}={\overline {{\bar {u}}_{i}{\bar {u}}_{j}}}-{\bar {u}}_{i}{\bar {u}}_{j}}" loading="lazy">, the Leonard tensor, represents interactions among large scales, <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle R_{ij}={\overline {u_{i}^{\prime }u_{j}^{\prime }}}}">

<msub>

</mrow>

</msub>

<mover>

<mrow>

</mrow>

</mrow>

</msubsup>

</mrow>

</mrow>

</msubsup>

</mrow>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle R_{ij}={\overline {u_{i}^{\prime }u_{j}^{\prime }}}}</annotation>

</semantics>

</math><img src="./43ae4844cc5dc88e47ef8b2d4b7eb253b56bbbf1.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.338ex; width:10.823ex; height:4.009ex;" alt="{\displaystyle R_{ij}={\overline {u_{i}^{\prime }u_{j}^{\prime }}}}" loading="lazy">, the Reynolds stress-like term, represents interactions among the sub-filter scales (SFS), and <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C_{ij}={\overline {{\bar {u}}_{i}u_{j}^{\prime }}}+{\overline {{\bar {u}}_{j}u_{i}^{\prime }}}}">

<msub>

</mrow>

</msub>

<mover>

<mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mrow>

</mrow>

</msubsup>

</mrow>

</mover>

</mrow>

<mover>

<mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mrow>

</mrow>

</msubsup>

</mrow>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C_{ij}={\overline {{\bar {u}}_{i}u_{j}^{\prime }}}+{\overline {{\bar {u}}_{j}u_{i}^{\prime }}}}</annotation>

</semantics>

</math><img src="./7124a8fdeecf6e4ed905046eb6167af9bca043c1.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.338ex; width:18.045ex; height:4.009ex;" alt="{\displaystyle C_{ij}={\overline {{\bar {u}}_{i}u_{j}^{\prime }}}+{\overline {{\bar {u}}_{j}u_{i}^{\prime }}}}" loading="lazy">, the Clark tensor,<a href="#cite_note-Clark-12">[12]</a> represents cross-scale interactions between large and small scales.<a href="#cite_note-Leonard_1974-11">[11]</a> Modeling the unclosed term <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \tau _{ij}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \tau _{ij}}</annotation>

</semantics>

</math><img src="./97d0190e7517254a655e8033badbc01c8590381c.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:2.493ex; height:2.343ex;" alt="{\displaystyle \tau _{ij}}" loading="lazy"> is the task of sub-grid scale (SGS) models. This is made challenging by the fact that the subgrid stress tensor <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \tau _{ij}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \tau _{ij}}</annotation>

</semantics>

</math><img src="./97d0190e7517254a655e8033badbc01c8590381c.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:2.493ex; height:2.343ex;" alt="{\displaystyle \tau _{ij}}" loading="lazy"> must account for interactions among all scales, including filtered scales with unfiltered scales.

The filtered governing equation for a passive scalar <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \phi }">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \phi }</annotation>

</semantics>

</math><img src="./72b1f30316670aee6270a28334bdf4f5072cdde4.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:1.385ex; height:2.509ex;" alt="{\displaystyle \phi }" loading="lazy">, such as mixture fraction or temperature, can be written as

<mfrac>

<mrow>

<mover>

</mover>

</mrow>

<mrow>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mrow>

<mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

<mover>

</mover>

</mrow>

</mrow>

<mfrac>

<mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\frac {\partial {\overline {\phi }}}{\partial t}}+{\frac {\partial }{\partial x_{j}}}\left({\overline {u}}_{j}{\overline {\phi }}\right)={\frac {\partial {\overline {J_{\phi }}}}{\partial x_{j}}}+{\frac {\partial q_{j}}{\partial x_{j}}}}</annotation>

</semantics>

</math><img src="./69bc42c7d1d7af731f284dc75351f4c0b53c9eaa.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.671ex; width:33.01ex; height:7.176ex;" alt="{\displaystyle {\frac {\partial {\overline {\phi }}}{\partial t}}+{\frac {\partial }{\partial x_{j}}}\left({\overline {u}}_{j}{\overline {\phi }}\right)={\frac {\partial {\overline {J_{\phi }}}}{\partial x_{j}}}+{\frac {\partial q_{j}}{\partial x_{j}}}}" loading="lazy"></dd></dl>

where <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle J_{\phi }}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle J_{\phi }}</annotation>

</semantics>

</math><img src="./84e8e06bcdc5ac75338af97d1a0c1df3f793fa73.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:2.502ex; height:2.843ex;" alt="{\displaystyle J_{\phi }}" loading="lazy"> is the diffusive flux of <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \phi }">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \phi }</annotation>

</semantics>

</math><img src="./72b1f30316670aee6270a28334bdf4f5072cdde4.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:1.385ex; height:2.509ex;" alt="{\displaystyle \phi }" loading="lazy">, and <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle q_{j}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle q_{j}}</annotation>

</semantics>

</math><img src="./e0d567ac2d170501680d2efa4c1d71d6a8569ef1.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:1.947ex; height:2.343ex;" alt="{\displaystyle q_{j}}" loading="lazy"> is the sub-filter flux for the scalar <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \phi }">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \phi }</annotation>

</semantics>

</math><img src="./72b1f30316670aee6270a28334bdf4f5072cdde4.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:1.385ex; height:2.509ex;" alt="{\displaystyle \phi }" loading="lazy">. The filtered diffusive flux <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\overline {J_{\phi }}}}">

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\overline {J_{\phi }}}}</annotation>

</semantics>

</math><img src="./4c11a4f94528c87e702af7dec29f57aa61c89f15.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:2.617ex; height:3.676ex;" alt="{\displaystyle {\overline {J_{\phi }}}}" loading="lazy"> is unclosed, unless a particular form is assumed for it, such as a gradient diffusion model <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle J_{\phi }=D_{\phi }{\frac {\partial \phi }{\partial x_{i}}}}">

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

<mfrac>

<mrow>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle J_{\phi }=D_{\phi }{\frac {\partial \phi }{\partial x_{i}}}}</annotation>

</semantics>

</math><img src="./0c5038055fb05492b05ef21734b1b2b6ec2bacff.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.338ex; width:13.02ex; height:6.009ex;" alt="{\displaystyle J_{\phi }=D_{\phi }{\frac {\partial \phi }{\partial x_{i}}}}" loading="lazy">. <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle q_{j}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle q_{j}}</annotation>

</semantics>

</math><img src="./e0d567ac2d170501680d2efa4c1d71d6a8569ef1.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:1.947ex; height:2.343ex;" alt="{\displaystyle q_{j}}" loading="lazy"> is defined analogously to <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \tau _{ij}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \tau _{ij}}</annotation>

</semantics>

</math><img src="./97d0190e7517254a655e8033badbc01c8590381c.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:2.493ex; height:2.343ex;" alt="{\displaystyle \tau _{ij}}" loading="lazy">,

<msub>

</mrow>

</msub>

<mover>

</mover>

</mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

<mover>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle q_{j}={\bar {\phi }}{\overline {u}}_{j}-{\overline {\phi u_{j}}}}</annotation>

</semantics>

</math><img src="./c07d657f60e7f9c15062fcdf88574e2cc4551c33.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:15.447ex; height:3.676ex;" alt="{\displaystyle q_{j}={\bar {\phi }}{\overline {u}}_{j}-{\overline {\phi u_{j}}}}" loading="lazy"></dd></dl>

and can similarly be split up into contributions from interactions between various scales. This sub-filter flux also requires a sub-filter model.

<div class="mw-heading mw-heading4"><h4 id="Derivation">Derivation</h4></div>

Using <a href="Einstein_notation" title="Einstein notation">Einstein notation</a>, the Navier–Stokes equations for an incompressible fluid in Cartesian coordinates are

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\frac {\partial u_{i}}{\partial x_{i}}}=0}</annotation>

</semantics>

</math><img src="./027a7e25a646c38061e3a3fbec1af683041abd81.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.338ex; width:8.544ex; height:5.843ex;" alt="{\displaystyle {\frac {\partial u_{i}}{\partial x_{i}}}=0}" loading="lazy"></dd>

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

<mrow>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

</mfrac>

</mrow>

<mfrac>

<mrow>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<msup>

</mrow>

</msup>

<msub>

</mrow>

</msub>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\frac {\partial u_{i}}{\partial t}}+{\frac {\partial u_{i}u_{j}}{\partial x_{j}}}=-{\frac {1}{\rho }}{\frac {\partial p}{\partial x_{i}}}+\nu {\frac {\partial ^{2}u_{i}}{\partial x_{j}\partial x_{j}}}.}</annotation>

</semantics>

</math><img src="./7dba4187a562b625c9e1484f5b2aa8b59181fd9b.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.671ex; width:37.546ex; height:6.509ex;" alt="{\displaystyle {\frac {\partial u_{i}}{\partial t}}+{\frac {\partial u_{i}u_{j}}{\partial x_{j}}}=-{\frac {1}{\rho }}{\frac {\partial p}{\partial x_{i}}}+\nu {\frac {\partial ^{2}u_{i}}{\partial x_{j}\partial x_{j}}}.}" loading="lazy"></dd></dl>

Filtering the momentum equation results in

<mover>

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

<mrow>

</mrow>

</mfrac>

</mover>

</mrow>

<mover>

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mover>

</mrow>

<mover>

<mrow>

<mfrac>

</mfrac>

</mrow>

<mfrac>

<mrow>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

</mover>

</mrow>

<mover>

<mrow>

<mfrac>

<mrow>

<msup>

</mrow>

</msup>

<msub>

</mrow>

</msub>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\overline {\frac {\partial u_{i}}{\partial t}}}+{\overline {\frac {\partial u_{i}u_{j}}{\partial x_{j}}}}=-{\overline {{\frac {1}{\rho }}{\frac {\partial p}{\partial x_{i}}}}}+{\overline {\nu {\frac {\partial ^{2}u_{i}}{\partial x_{j}\partial x_{j}}}}}.}</annotation>

</semantics>

</math><img src="./fff82183eb3e3c5ed873df4d83541f421c9a045c.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.671ex; width:38.005ex; height:7.176ex;" alt="{\displaystyle {\overline {\frac {\partial u_{i}}{\partial t}}}+{\overline {\frac {\partial u_{i}u_{j}}{\partial x_{j}}}}=-{\overline {{\frac {1}{\rho }}{\frac {\partial p}{\partial x_{i}}}}}+{\overline {\nu {\frac {\partial ^{2}u_{i}}{\partial x_{j}\partial x_{j}}}}}.}" loading="lazy"></dd></dl>

If we assume that filtering and differentiation commute, then

<mfrac>

<mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mrow>

</mrow>

</mfrac>

</mrow>

<mover>

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mover>

</mrow>

<mfrac>

</mfrac>

</mrow>

<mfrac>

<mrow>

<mover>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<msup>

</mrow>

</msup>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\frac {\partial {\bar {u_{i}}}}{\partial t}}+{\overline {\frac {\partial u_{i}u_{j}}{\partial x_{j}}}}=-{\frac {1}{\rho }}{\frac {\partial {\bar {p}}}{\partial x_{i}}}+\nu {\frac {\partial ^{2}{\bar {u_{i}}}}{\partial x_{j}\partial x_{j}}}.}</annotation>

</semantics>

</math><img src="./47cacaf1095658cedb266e51af017782df58fa8e.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.671ex; width:37.661ex; height:7.176ex;" alt="{\displaystyle {\frac {\partial {\bar {u_{i}}}}{\partial t}}+{\overline {\frac {\partial u_{i}u_{j}}{\partial x_{j}}}}=-{\frac {1}{\rho }}{\frac {\partial {\bar {p}}}{\partial x_{i}}}+\nu {\frac {\partial ^{2}{\bar {u_{i}}}}{\partial x_{j}\partial x_{j}}}.}" loading="lazy"></dd></dl>

This equation models the changes in time of the filtered variables <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\bar {u_{i}}}}">

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\bar {u_{i}}}}</annotation>

</semantics>

</math><img src="./7702717043d411c5fdd38599f7ff20751df794bc.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:2.129ex; height:2.343ex;" alt="{\displaystyle {\bar {u_{i}}}}" loading="lazy">. Since the unfiltered variables <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle u_{i}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle u_{i}}</annotation>

</semantics>

</math><img src="./14f13cb025ff2e136dcbd2fc81ddf965b728e3d7.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:2.129ex; height:2.009ex;" alt="{\displaystyle u_{i}}" loading="lazy"> are not known, it is impossible to directly calculate <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\overline {\frac {\partial u_{i}u_{j}}{\partial x_{j}}}}}">

<mover>

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\overline {\frac {\partial u_{i}u_{j}}{\partial x_{j}}}}}</annotation>

</semantics>

</math><img src="./3704aa5eed681b197a5c16e950e9e0bc13aba739.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.671ex; width:6.638ex; height:7.176ex;" alt="{\displaystyle {\overline {\frac {\partial u_{i}u_{j}}{\partial x_{j}}}}}" loading="lazy">. However, the quantity <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\frac {\partial {\bar {u_{i}}}{\bar {u_{j}}}}{\partial x_{j}}}}">

<mfrac>

<mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\frac {\partial {\bar {u_{i}}}{\bar {u_{j}}}}{\partial x_{j}}}}</annotation>

</semantics>

</math><img src="./40710a80b93f98966872a01e6b847bee6f697ef6.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.671ex; width:6.523ex; height:6.509ex;" alt="{\displaystyle {\frac {\partial {\bar {u_{i}}}{\bar {u_{j}}}}{\partial x_{j}}}}" loading="lazy"> is known. A substitution is made:

<mfrac>

<mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mrow>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

</mfrac>

</mrow>

<mfrac>

<mrow>

<mover>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<msup>

</mrow>

</msup>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mrow>

<mrow>

<mover>

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mover>

</mrow>

<mfrac>

<mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\frac {\partial {\bar {u_{i}}}}{\partial t}}+{\frac {\partial {\bar {u_{i}}}{\bar {u_{j}}}}{\partial x_{j}}}=-{\frac {1}{\rho }}{\frac {\partial {\bar {p}}}{\partial x_{i}}}+\nu {\frac {\partial ^{2}{\bar {u_{i}}}}{\partial x_{j}\partial x_{j}}}-\left({\overline {\frac {\partial u_{i}u_{j}}{\partial x_{j}}}}-{\frac {\partial {\bar {u_{i}}}{\bar {u_{j}}}}{\partial x_{j}}}\right).}</annotation>

</semantics>

</math><img src="./6876ec5a46f72a129750ac51a5f888aa2a5d9e13.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -3.171ex; width:60.456ex; height:7.676ex;" alt="{\displaystyle {\frac {\partial {\bar {u_{i}}}}{\partial t}}+{\frac {\partial {\bar {u_{i}}}{\bar {u_{j}}}}{\partial x_{j}}}=-{\frac {1}{\rho }}{\frac {\partial {\bar {p}}}{\partial x_{i}}}+\nu {\frac {\partial ^{2}{\bar {u_{i}}}}{\partial x_{j}\partial x_{j}}}-\left({\overline {\frac {\partial u_{i}u_{j}}{\partial x_{j}}}}-{\frac {\partial {\bar {u_{i}}}{\bar {u_{j}}}}{\partial x_{j}}}\right).}" loading="lazy"></dd></dl>

Let <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \tau _{ij}={\overline {u_{i}u_{j}}}-{\bar {u}}_{i}{\bar {u}}_{j}}">

<msub>

</mrow>

</msub>

<mover>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

</mover>

</mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \tau _{ij}={\overline {u_{i}u_{j}}}-{\bar {u}}_{i}{\bar {u}}_{j}}</annotation>

</semantics>

</math><img src="./15e2e14ac54b637d683e874dc63e74d15dbaef51.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:17.285ex; height:3.009ex;" alt="{\displaystyle \tau _{ij}={\overline {u_{i}u_{j}}}-{\bar {u}}_{i}{\bar {u}}_{j}}" loading="lazy">. The resulting set of equations are the LES equations:

<mfrac>

<mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mrow>

</mrow>

</mfrac>

</mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mfrac>

<mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

</mfrac>

</mrow>

<mfrac>

<mrow>

<mover>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<msup>

</mrow>

</msup>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\frac {\partial {\bar {u_{i}}}}{\partial t}}+{\bar {u_{j}}}{\frac {\partial {\bar {u_{i}}}}{\partial x_{j}}}=-{\frac {1}{\rho }}{\frac {\partial {\bar {p}}}{\partial x_{i}}}+\nu {\frac {\partial ^{2}{\bar {u_{i}}}}{\partial x_{j}\partial x_{j}}}-{\frac {\partial \tau _{ij}}{\partial x_{j}}}.}</annotation>

</semantics>

</math><img src="./b5a78c791052a2a42af185622ab52c6a562cca0c.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.671ex; width:45.143ex; height:6.509ex;" alt="{\displaystyle {\frac {\partial {\bar {u_{i}}}}{\partial t}}+{\bar {u_{j}}}{\frac {\partial {\bar {u_{i}}}}{\partial x_{j}}}=-{\frac {1}{\rho }}{\frac {\partial {\bar {p}}}{\partial x_{i}}}+\nu {\frac {\partial ^{2}{\bar {u_{i}}}}{\partial x_{j}\partial x_{j}}}-{\frac {\partial \tau _{ij}}{\partial x_{j}}}.}" loading="lazy"></dd></dl>

<div class="mw-heading mw-heading3"><h3 id="Compressible_governing_equations">Compressible governing equations</h3></div>

For the governing equations of compressible flow, each equation, starting with the conservation of mass, is filtered. This gives:

<mfrac>

<mrow>

<mover>

</mover>

</mrow>

<mrow>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<mover>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\frac {\partial {\overline {\rho }}}{\partial t}}+{\frac {\partial {\overline {u_{i}\rho }}}{\partial x_{i}}}=0}</annotation>

</semantics>

</math><img src="./30f735aea67cdc8e70face6a5960fb87dc04eb01.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.338ex; width:16.173ex; height:6.176ex;" alt="{\displaystyle {\frac {\partial {\overline {\rho }}}{\partial t}}+{\frac {\partial {\overline {u_{i}\rho }}}{\partial x_{i}}}=0}" loading="lazy"></dd></dl>

which results in an additional sub-filter term. However, it is desirable to avoid having to model the sub-filter scales of the mass conservation equation. For this reason, Favre<a href="#cite_note-Favre_1983-13">[13]</a> proposed a density-weighted filtering operation, called Favre filtering, defined for an arbitrary quantity <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \phi }">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \phi }</annotation>

</semantics>

</math><img src="./72b1f30316670aee6270a28334bdf4f5072cdde4.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:1.385ex; height:2.509ex;" alt="{\displaystyle \phi }" loading="lazy"> as:

<mover>

</mover>

</mrow>

<mfrac>

<mover>

<mrow>

</mrow>

</mover>

<mover>

</mover>

</mfrac>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\tilde {\phi }}={\frac {\overline {\rho \phi }}{\overline {\rho }}}}</annotation>

</semantics>

</math><img src="./a1467eb42e8f4321bc3359e842471fda4bde1e1c.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.671ex; width:8.104ex; height:7.009ex;" alt="{\displaystyle {\tilde {\phi }}={\frac {\overline {\rho \phi }}{\overline {\rho }}}}" loading="lazy"></dd></dl>

which, in the limit of incompressibility, becomes the normal filtering operation. This makes the conservation of mass equation:

<mfrac>

<mrow>

<mover>

</mover>

</mrow>

<mrow>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<mover>

</mover>

</mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\frac {\partial {\overline {\rho }}}{\partial t}}+{\frac {\partial {\overline {\rho }}{\tilde {u_{i}}}}{\partial x_{i}}}=0.}</annotation>

</semantics>

</math><img src="./5631754c633fa0665ed725d1188cba548db50650.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.338ex; width:16.82ex; height:6.176ex;" alt="{\displaystyle {\frac {\partial {\overline {\rho }}}{\partial t}}+{\frac {\partial {\overline {\rho }}{\tilde {u_{i}}}}{\partial x_{i}}}=0.}" loading="lazy"></dd></dl>

This concept can then be extended to write the Favre-filtered momentum equation for compressible flow. Following Vreman:<a href="#cite_note-Vreman_1995-14">[14]</a>

<mfrac>

<mrow>

<mover>

</mover>

</mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mrow>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<mover>

</mover>

</mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<mover>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<mover>

</mover>

</mrow>

</mrow>

</mrow>

</msubsup>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mrow>

<mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mrow>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\frac {\partial {\overline {\rho }}{\tilde {u_{i}}}}{\partial t}}+{\frac {\partial {\overline {\rho }}{\tilde {u_{i}}}{\tilde {u_{j}}}}{\partial x_{j}}}+{\frac {\partial {\overline {p}}}{\partial x_{i}}}-{\frac {\partial {\tilde {\sigma }}_{ij}}{\partial x_{j}}}=-{\frac {\partial {\overline {\rho }}\tau _{ij}^{r}}{\partial x_{j}}}+{\frac {\partial }{\partial x_{j}}}\left({\overline {\sigma }}_{ij}-{\tilde {\sigma }}_{ij}\right)}</annotation>

</semantics>

</math><img src="./9c5c64e44e4470dc0fc784a307e0dcc39abf6e34.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.671ex; width:60.078ex; height:7.009ex;" alt="{\displaystyle {\frac {\partial {\overline {\rho }}{\tilde {u_{i}}}}{\partial t}}+{\frac {\partial {\overline {\rho }}{\tilde {u_{i}}}{\tilde {u_{j}}}}{\partial x_{j}}}+{\frac {\partial {\overline {p}}}{\partial x_{i}}}-{\frac {\partial {\tilde {\sigma }}_{ij}}{\partial x_{j}}}=-{\frac {\partial {\overline {\rho }}\tau _{ij}^{r}}{\partial x_{j}}}+{\frac {\partial }{\partial x_{j}}}\left({\overline {\sigma }}_{ij}-{\tilde {\sigma }}_{ij}\right)}" loading="lazy"></dd></dl>

where <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \sigma _{ij}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \sigma _{ij}}</annotation>

</semantics>

</math><img src="./43acbf52cc4d4f83f187ceaa49f045114b71772e.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:2.804ex; height:2.343ex;" alt="{\displaystyle \sigma _{ij}}" loading="lazy"> is the <a href="Shear_stress" title="Shear stress">shear stress</a> tensor, given for a <a href="Newtonian_fluid" title="Newtonian fluid">Newtonian fluid</a> by:

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

<mfrac>

</mfrac>

</mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \sigma _{ij}=2\mu (T)S_{ij}-{\frac {2}{3}}\mu (T)\delta _{ij}S_{kk}}</annotation>

</semantics>

</math><img src="./82b7f201ec91abf62f80e8f373da85f01fed0a9a.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.838ex; width:30.381ex; height:5.176ex;" alt="{\displaystyle \sigma _{ij}=2\mu (T)S_{ij}-{\frac {2}{3}}\mu (T)\delta _{ij}S_{kk}}" loading="lazy"></dd></dl>

and the term <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\frac {\partial }{\partial x_{j}}}\left({\overline {\sigma }}_{ij}-{\tilde {\sigma }}_{ij}\right)}">

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mrow>

<mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mrow>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\frac {\partial }{\partial x_{j}}}\left({\overline {\sigma }}_{ij}-{\tilde {\sigma }}_{ij}\right)}</annotation>

</semantics>

</math><img src="./2935d447f64d571ddc1356f4916ef1348653cbaf.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.671ex; width:15.161ex; height:6.176ex;" alt="{\displaystyle {\frac {\partial }{\partial x_{j}}}\left({\overline {\sigma }}_{ij}-{\tilde {\sigma }}_{ij}\right)}" loading="lazy"> represents a sub-filter viscous contribution from evaluating the viscosity <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \mu (T)}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \mu (T)}</annotation>

</semantics>

</math><img src="./a58da74e9341eac12fc425f0eacd0f83a5183711.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:4.847ex; height:2.843ex;" alt="{\displaystyle \mu (T)}" loading="lazy"> using the Favre-filtered temperature <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\tilde {T}}}">

<mover>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\tilde {T}}}</annotation>

</semantics>

</math><img src="./4827ccc5dba8bc62d9ea7a3d129cd3915582fa89.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.774ex; height:2.676ex;" alt="{\displaystyle {\tilde {T}}}" loading="lazy">. The subgrid stress tensor for the Favre-filtered momentum field is given by

</mrow>

</mrow>

</msubsup>

<mover>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

</mover>

</mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \tau _{ij}^{r}={\widetilde {u_{i}\cdot u_{j}}}-{\tilde {u_{i}}}{\tilde {u_{j}}}}</annotation>

</semantics>

</math><img src="./d6068fe665a8172b2040a284ba1d1ae073bd2ae1.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.338ex; width:18.849ex; height:3.509ex;" alt="{\displaystyle \tau _{ij}^{r}={\widetilde {u_{i}\cdot u_{j}}}-{\tilde {u_{i}}}{\tilde {u_{j}}}}" loading="lazy"></dd></dl>

By analogy, the Leonard decomposition may also be written for the residual stress tensor for a filtered triple product <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\overline {\rho \phi \psi }}}">

<mover>

<mrow>

</mrow>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\overline {\rho \phi \psi }}}</annotation>

</semantics>

</math><img src="./028af53c3e666a656379d73c8129d8d31b56c564.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:4.215ex; height:3.509ex;" alt="{\displaystyle {\overline {\rho \phi \psi }}}" loading="lazy">. The triple product can be rewritten using the Favre filtering operator as <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\overline {\rho }}{\widetilde {\phi \psi }}}">

<mover>

</mover>

</mrow>

<mover>

<mrow>

</mrow>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\overline {\rho }}{\widetilde {\phi \psi }}}</annotation>

</semantics>

</math><img src="./1075ef34cdb48ac95a0ad7da03ef89aae0e73824.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:4.215ex; height:3.343ex;" alt="{\displaystyle {\overline {\rho }}{\widetilde {\phi \psi }}}" loading="lazy">, which is an unclosed term (it requires knowledge of the fields <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \phi }">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \phi }</annotation>

</semantics>

</math><img src="./72b1f30316670aee6270a28334bdf4f5072cdde4.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:1.385ex; height:2.509ex;" alt="{\displaystyle \phi }" loading="lazy"> and <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \psi }">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \psi }</annotation>

</semantics>

</math><img src="./45e5789e5d9c8f7c79744f43ecaaf8ba42a8553a.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:1.513ex; height:2.509ex;" alt="{\displaystyle \psi }" loading="lazy">, when only the fields <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\tilde {\phi }}}">

<mover>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\tilde {\phi }}}</annotation>

</semantics>

</math><img src="./d9d195019ecfbb56473440008ed4b41a7ee46fe0.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:1.467ex; height:3.009ex;" alt="{\displaystyle {\tilde {\phi }}}" loading="lazy"> and <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\tilde {\psi }}}">

<mover>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\tilde {\psi }}}</annotation>

</semantics>

</math><img src="./a4ac6d8bfc5b19aea7fd59edc84dca2a1fa3a86e.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:1.596ex; height:3.009ex;" alt="{\displaystyle {\tilde {\psi }}}" loading="lazy"> are known). It can be broken up in a manner analogous to <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\overline {u_{i}u_{j}}}}">

<mover>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\overline {u_{i}u_{j}}}}</annotation>

</semantics>

</math><img src="./a850035c59b486247399b36c5ef57af0a4c88a37.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:4.484ex; height:3.009ex;" alt="{\displaystyle {\overline {u_{i}u_{j}}}}" loading="lazy"> above, which results in a sub-filter stress tensor <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\overline {\rho }}\left({\widetilde {\phi \psi }}-{\tilde {\phi }}{\tilde {\psi }}\right)}">

<mover>

</mover>

</mrow>

<mrow>

<mrow>

<mover>

<mrow>

</mrow>

</mover>

</mrow>

<mover>

</mover>

</mrow>

<mover>

</mover>

</mrow>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\overline {\rho }}\left({\widetilde {\phi \psi }}-{\tilde {\phi }}{\tilde {\psi }}\right)}</annotation>

</semantics>

</math><img src="./79df9a0671f6c56ec605727281053d68aade9cde.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.838ex; width:13.282ex; height:4.843ex;" alt="{\displaystyle {\overline {\rho }}\left({\widetilde {\phi \psi }}-{\tilde {\phi }}{\tilde {\psi }}\right)}" loading="lazy">. This sub-filter term can be split up into contributions from three types of interactions: the Leondard tensor <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle L_{ij}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle L_{ij}}</annotation>

</semantics>

</math><img src="./a310c7a19bc362591bb6809bb95b7e67baa59965.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:3.06ex; height:2.843ex;" alt="{\displaystyle L_{ij}}" loading="lazy">, representing interactions among resolved scales; the Clark tensor <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C_{ij}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C_{ij}}</annotation>

</semantics>

</math><img src="./70b3f8ea18a9b9a0dbe09dae212a623820b217d7.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:3.139ex; height:2.843ex;" alt="{\displaystyle C_{ij}}" loading="lazy">, representing interactions between resolved and unresolved scales; and the Reynolds tensor <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle R_{ij}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle R_{ij}}</annotation>

</semantics>

</math><img src="./50382df1fd3bebf0baf22e2eab56a7a471d6621a.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:3.241ex; height:2.843ex;" alt="{\displaystyle R_{ij}}" loading="lazy">, which represents interactions among unresolved scales.<a href="#cite_note-Sagaut_2009-15">[15]</a>

<div class="mw-heading mw-heading3"><h3 id="Filtered_kinetic_energy_equation">Filtered kinetic energy equation</h3></div>

In addition to the filtered mass and momentum equations, filtering the kinetic energy equation can provide additional insight. The kinetic energy field can be filtered to yield the total filtered kinetic energy:

<mover>

</mover>

</mrow>

<mfrac>

</mfrac>

</mrow>

<mover>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\overline {E}}={\frac {1}{2}}{\overline {u_{i}u_{i}}}}</annotation>

</semantics>

</math><img src="./e3b96f700f9f3d29c02b765744568d6ffc01a1e3.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.838ex; width:11.401ex; height:5.176ex;" alt="{\displaystyle {\overline {E}}={\frac {1}{2}}{\overline {u_{i}u_{i}}}}" loading="lazy"></dd></dl>

and the total filtered kinetic energy can be decomposed into two terms: the kinetic energy of the filtered velocity field <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle E_{f}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle E_{f}}</annotation>

</semantics>

</math><img src="./12ce471791c1904f7401133fb39c86f5a088f205.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:2.852ex; height:2.843ex;" alt="{\displaystyle E_{f}}" loading="lazy">,

<msub>

</mrow>

</msub>

<mfrac>

</mfrac>

</mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle E_{f}={\frac {1}{2}}{\overline {u_{i}}}\,{\overline {u_{i}}}}</annotation>

</semantics>

</math><img src="./20f605b516798ae48cf9e7662b57deacb7b2f226.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.838ex; width:12.824ex; height:5.176ex;" alt="{\displaystyle E_{f}={\frac {1}{2}}{\overline {u_{i}}}\,{\overline {u_{i}}}}" loading="lazy"></dd></dl>

and the residual kinetic energy <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle k_{r}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle k_{r}}</annotation>

</semantics>

</math><img src="./58d7e5bb94a98612867e7ad219e44b5499bb7ebe.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:2.185ex; height:2.509ex;" alt="{\displaystyle k_{r}}" loading="lazy">,

<msub>

</mrow>

</msub>

<mfrac>

</mfrac>

</mrow>

<mover>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

</mover>

</mrow>

<mfrac>

</mfrac>

</mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mfrac>

</mfrac>

</mrow>

</mrow>

</mrow>

</msubsup>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle k_{r}={\frac {1}{2}}{\overline {u_{i}u_{i}}}-{\frac {1}{2}}{\overline {u_{i}}}\,{\overline {u_{i}}}={\frac {1}{2}}\tau _{ii}^{r}}</annotation>

</semantics>

</math><img src="./5a76582fc30be63ade0bf239a738505d5d8d610a.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.838ex; width:28.851ex; height:5.176ex;" alt="{\displaystyle k_{r}={\frac {1}{2}}{\overline {u_{i}u_{i}}}-{\frac {1}{2}}{\overline {u_{i}}}\,{\overline {u_{i}}}={\frac {1}{2}}\tau _{ii}^{r}}" loading="lazy"></dd></dl>

such that <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\overline {E}}=E_{f}+k_{r}}">

<mover>

</mover>

</mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\overline {E}}=E_{f}+k_{r}}</annotation>

</semantics>

</math><img src="./46d6911866f2fa98cc2184c9872ded6b64d0d6c2.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:12.905ex; height:3.676ex;" alt="{\displaystyle {\overline {E}}=E_{f}+k_{r}}" loading="lazy">.

The conservation equation for <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle E_{f}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle E_{f}}</annotation>

</semantics>

</math><img src="./12ce471791c1904f7401133fb39c86f5a088f205.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:2.852ex; height:2.843ex;" alt="{\displaystyle E_{f}}" loading="lazy"> can be obtained by multiplying the filtered momentum transport equation by <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\overline {u_{i}}}}">

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\overline {u_{i}}}}</annotation>

</semantics>

</math><img src="./2be16dfddf2b18a82a6b3f026e55f5e4ee9a91e8.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:2.244ex; height:2.676ex;" alt="{\displaystyle {\overline {u_{i}}}}" loading="lazy"> to yield:

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

<mrow>

</mrow>

</mfrac>

</mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

</mfrac>

</mrow>

<mfrac>

<mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mover>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

</mrow>

</mrow>

</msubsup>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\frac {\partial E_{f}}{\partial t}}+{\overline {u_{j}}}{\frac {\partial E_{f}}{\partial x_{j}}}+{\frac {1}{\rho }}{\frac {\partial {\overline {u_{i}}}{\bar {p}}}{\partial x_{i}}}+{\frac {\partial {\overline {u_{i}}}\tau _{ij}^{r}}{\partial x_{j}}}-2\nu {\frac {\partial {\overline {u_{i}}}{\bar {S_{ij}}}}{\partial x_{j}}}=-\epsilon _{f}-\Pi }</annotation>

</semantics>

</math><img src="./8e17e72549ac832f6965a527b7a07e81e5a7c7fb.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.671ex; width:59.681ex; height:7.009ex;" alt="{\displaystyle {\frac {\partial E_{f}}{\partial t}}+{\overline {u_{j}}}{\frac {\partial E_{f}}{\partial x_{j}}}+{\frac {1}{\rho }}{\frac {\partial {\overline {u_{i}}}{\bar {p}}}{\partial x_{i}}}+{\frac {\partial {\overline {u_{i}}}\tau _{ij}^{r}}{\partial x_{j}}}-2\nu {\frac {\partial {\overline {u_{i}}}{\bar {S_{ij}}}}{\partial x_{j}}}=-\epsilon _{f}-\Pi }" loading="lazy"></dd></dl>

where <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \epsilon _{f}=2\nu {\bar {S_{ij}}}{\bar {S_{ij}}}}">

<msub>

</mrow>

</msub>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \epsilon _{f}=2\nu {\bar {S_{ij}}}{\bar {S_{ij}}}}</annotation>

</semantics>

</math><img src="./67cd85f4752f3e50d10fa9bfb050fea06eb59bb0.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:13.378ex; height:3.343ex;" alt="{\displaystyle \epsilon _{f}=2\nu {\bar {S_{ij}}}{\bar {S_{ij}}}}" loading="lazy"> is the dissipation of kinetic energy of the filtered velocity field by viscous stress, and <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \Pi =-\tau _{ij}^{r}{\bar {S_{ij}}}}">

</mrow>

</mrow>

</msubsup>

<mover>

<msub>

</mrow>

</msub>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \Pi =-\tau _{ij}^{r}{\bar {S_{ij}}}}</annotation>

</semantics>

</math><img src="./616881bdb2194427c1cb1b31074695842932cfe2.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.338ex; width:12.045ex; height:3.676ex;" alt="{\displaystyle \Pi =-\tau _{ij}^{r}{\bar {S_{ij}}}}" loading="lazy"> represents the sub-filter scale (SFS) dissipation of kinetic energy.

The terms on the left-hand side represent transport, and the terms on the right-hand side are sink terms that dissipate kinetic energy.<a href="#cite_note-Pope_2000-9">[9]</a>

The <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \Pi }">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \Pi }</annotation>

</semantics>

</math><img src="./eed3e3db6cc2028a183af948212ed2551d25c954.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.743ex; height:2.176ex;" alt="{\displaystyle \Pi }" loading="lazy"> SFS dissipation term is of particular interest, since it represents the transfer of energy from large resolved scales to small unresolved scales. On average, <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \Pi }">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \Pi }</annotation>

</semantics>

</math><img src="./eed3e3db6cc2028a183af948212ed2551d25c954.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.743ex; height:2.176ex;" alt="{\displaystyle \Pi }" loading="lazy"> transfers energy from large to small scales. However, instantaneously <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \Pi }">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \Pi }</annotation>

</semantics>

</math><img src="./eed3e3db6cc2028a183af948212ed2551d25c954.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.743ex; height:2.176ex;" alt="{\displaystyle \Pi }" loading="lazy"> can be positive or negative, meaning it can also act as a source term for <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle E_{f}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle E_{f}}</annotation>

</semantics>

</math><img src="./12ce471791c1904f7401133fb39c86f5a088f205.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:2.852ex; height:2.843ex;" alt="{\displaystyle E_{f}}" loading="lazy">, the kinetic energy of the filtered velocity field. The transfer of energy from unresolved to resolved scales is called backscatter (and likewise the transfer of energy from resolved to unresolved scales is called forward-scatter).<a href="#cite_note-Piomelli_1991-16">[16]</a>

<div class="mw-heading mw-heading2"><h2 id="Numerical_methods_for_LES">Numerical methods for LES</h2></div>

Large eddy simulation involves the solution to the discrete filtered governing equations using <a href="Computational_fluid_dynamics" title="Computational fluid dynamics">computational fluid dynamics</a>. LES resolves scales from the domain size <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle L}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle L}</annotation>

</semantics>

</math><img src="./103168b86f781fe6e9a4a87b8ea1cebe0ad4ede8.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.583ex; height:2.176ex;" alt="{\displaystyle L}" loading="lazy"> down to the filter size <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \Delta }">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \Delta }</annotation>

</semantics>

</math><img src="./32769037c408874e1890f77554c65f39c523ebe2.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.936ex; height:2.176ex;" alt="{\displaystyle \Delta }" loading="lazy">, and as such a substantial portion of high wave number turbulent fluctuations must be resolved. This requires either <a href="High-resolution_scheme" title="High-resolution scheme">high-order numerical schemes</a>, or fine grid resolution if low-order numerical schemes are used. Chapter 13 of Pope<a href="#cite_note-Pope_2000-9">[9]</a> addresses the question of how fine a grid resolution <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \Delta x}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \Delta x}</annotation>

</semantics>

</math><img src="./f3890eb866b6258d7a304fc34c70ee3fb3a81a70.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:3.266ex; height:2.176ex;" alt="{\displaystyle \Delta x}" loading="lazy"> is needed to resolve a filtered velocity field <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\overline {u}}({\boldsymbol {x}})}">

<mover>

</mover>

</mrow>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\overline {u}}({\boldsymbol {x}})}</annotation>

</semantics>

</math><img src="./e85175ef82d4db557c9f6ced918e4d6cf10575ae.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:4.786ex; height:2.843ex;" alt="{\displaystyle {\overline {u}}({\boldsymbol {x}})}" loading="lazy">. Ghosal<a href="#cite_note-Ghosal_1996-17">[17]</a> found that for low-order discretization schemes, such as those used in finite volume methods, the truncation error can be the same order as the subfilter scale contributions, unless the filter width <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \Delta }">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \Delta }</annotation>

</semantics>

</math><img src="./32769037c408874e1890f77554c65f39c523ebe2.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.936ex; height:2.176ex;" alt="{\displaystyle \Delta }" loading="lazy"> is considerably larger than the grid spacing <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \Delta x}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \Delta x}</annotation>

</semantics>

</math><img src="./f3890eb866b6258d7a304fc34c70ee3fb3a81a70.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:3.266ex; height:2.176ex;" alt="{\displaystyle \Delta x}" loading="lazy">. While even-order schemes have truncation error, they are non-dissipative,<a href="#cite_note-Leveque_1992-18">[18]</a> and because subfilter scale models are dissipative, even-order schemes will not affect the subfilter scale model contributions as strongly as dissipative schemes.

<div class="mw-heading mw-heading3"><h3 id="Filter_implementation">Filter implementation</h3></div>

The filtering operation in large eddy simulation can be implicit or explicit. Implicit filtering recognizes that the subfilter scale model will dissipate in the same manner as many numerical schemes. In this way, the grid, or the numerical discretization scheme, can be assumed to be the LES low-pass filter. While this takes full advantage of the grid resolution, and eliminates the computational cost of calculating a subfilter scale model term, it is difficult to determine the shape of the LES filter that is associated with some numerical issues. Additionally, truncation error can also become an issue.<a href="#cite_note-Grinstein_2007-19">[19]</a>

In explicit filtering, an <a href="Filter_(large_eddy_simulation)" title="Filter (large eddy simulation)">LES filter</a> is applied to the discretized Navier–Stokes equations, providing a well-defined filter shape and reducing the truncation error. However, explicit filtering requires a finer grid than implicit filtering, and the computational cost increases with <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle (\Delta x)^{4}}">

<msup>

</mrow>

</msup>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle (\Delta x)^{4}}</annotation>

</semantics>

</math><img src="./7f4a3bf4ed873c87a817a6b7f9e387b57faa67f7.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:6.129ex; height:3.176ex;" alt="{\displaystyle (\Delta x)^{4}}" loading="lazy">. Chapter 8 of Sagaut (2006) covers LES numerics in greater detail.<a href="#cite_note-Sagaut_2006-10">[10]</a>

<div class="mw-heading mw-heading2"><h2 id="Boundary_conditions_of_large_eddy_simulations">Boundary conditions of large eddy simulations</h2></div>

Inlet boundary conditions affect the accuracy of LES significantly, and the treatment of inlet conditions for LES is a complicated problem. Theoretically, a good boundary condition for LES should contain the following features:<a href="#cite_note-Li_P-20">[20]</a>

(1) providing accurate information of flow characteristics, i.e. velocity and turbulence;

(2) satisfying the Navier-Stokes equations and other physics;

(3) being easy to implement and adjust to different cases.

Currently, methods of generating inlet conditions for LES are broadly divided into two categories classified by Tabor et al.:<a href="#cite_note-21">[21]</a>

The first method for generating turbulent inlets is to synthesize them according to particular cases, such as Fourier techniques, principle orthogonal decomposition (POD) and vortex methods. The synthesis techniques attempt to construct turbulent field at inlets that have suitable turbulence-like properties and make it easy to specify parameters of the turbulence, such as turbulent kinetic energy and turbulent dissipation rate. In addition, inlet conditions generated by using random numbers are computationally inexpensive. However, one serious drawback exists in the method. The synthesized turbulence does not satisfy the physical structure of fluid flow governed by Navier-Stokes equations.<a href="#cite_note-Li_P-20">[20]</a>

The second method involves a separate and precursor calculation to generate a turbulent database which can be introduced into the main computation at the inlets. The database (sometimes named as ‘library’) can be generated in a number of ways, such as cyclic domains, pre-prepared library, and internal mapping. However, the method of generating turbulent inflow by precursor simulations requires large calculation capacity.

Researchers examining the application of various types of synthetic and precursor calculations have found that the more realistic the inlet turbulence, the more accurate LES predicts results.<a href="#cite_note-Li_P-20">[20]</a>

<div class="mw-heading mw-heading2"><h2 id="Modeling_unresolved_scales">Modeling unresolved scales</h2></div>

To discuss the modeling of unresolved scales, first the unresolved scales must be classified. They fall into two groups: resolved sub-filter scales (SFS), and sub-grid scales(SGS).

The resolved sub-filter scales represent the scales with wave numbers larger than the cutoff wave number <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle k_{c}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle k_{c}}</annotation>

</semantics>

</math><img src="./0d11e2ff9fda36c0dd027d4f15f3ce1607c037fd.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:2.155ex; height:2.509ex;" alt="{\displaystyle k_{c}}" loading="lazy">, but whose effects are dampened by the filter. Resolved sub-filter scales only exist when filters non-local in wave-space are used (such as a <a href="Filter_(large_eddy_simulation)#Box_filter" title="Filter (large eddy simulation)">box</a> or <a href="Filter_(large_eddy_simulation)#Gaussian_filter" title="Filter (large eddy simulation)">Gaussian</a> filter). These resolved sub-filter scales must be modeled using filter reconstruction.

Sub-grid scales are any scales that are smaller than the cutoff filter width <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \Delta }">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \Delta }</annotation>

</semantics>

</math><img src="./32769037c408874e1890f77554c65f39c523ebe2.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.936ex; height:2.176ex;" alt="{\displaystyle \Delta }" loading="lazy">. The form of the SGS model depends on the filter implementation. As mentioned in the <a href="#Numerical_Methods_for_LES">Numerical methods for LES</a> section, if implicit LES is considered, no SGS model is implemented and the numerical effects of the discretization are assumed to mimic the physics of the unresolved turbulent motions.

<div class="mw-heading mw-heading3"><h3 id="Sub-grid_scale_models">Sub-grid scale models</h3></div>

Without a universally valid description of turbulence, empirical information must be utilized when constructing and applying SGS models, supplemented with fundamental physical constraints such as <a href="Galilean_invariance" title="Galilean invariance">Galilean invariance</a><a href="#cite_note-Pope_2000-9">[9]</a>

.<a href="#cite_note-Meneveau_2010-22">[22]</a>

Two classes of SGS models exist; the first class is functional models and the second class is structural models. Some models may be categorized as both.

<div class="mw-heading mw-heading4"><h4 id="Functional_(eddy–viscosity)_models">Functional (eddy–viscosity) models</h4></div>

Functional models are simpler than structural models, focusing only on dissipating energy at a rate that is physically correct. These are based on an artificial eddy viscosity approach, where the effects of turbulence are lumped into a turbulent viscosity. The approach treats dissipation of kinetic energy at sub-grid scales as analogous to molecular diffusion. In this case, the deviatoric part of <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \tau _{ij}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \tau _{ij}}</annotation>

</semantics>

</math><img src="./97d0190e7517254a655e8033badbc01c8590381c.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:2.493ex; height:2.343ex;" alt="{\displaystyle \tau _{ij}}" loading="lazy"> is modeled as:

</mrow>

</mrow>

</msubsup>

<mfrac>

</mfrac>

</mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \tau _{ij}^{r}-{\frac {1}{3}}\tau _{kk}\delta _{ij}=-2\nu _{\mathrm {t} }{\bar {S}}_{ij}}</annotation>

</semantics>

</math><img src="./c39d01001589ef92ac8e6cf7bd11152f280faf47.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.838ex; width:23.942ex; height:5.176ex;" alt="{\displaystyle \tau _{ij}^{r}-{\frac {1}{3}}\tau _{kk}\delta _{ij}=-2\nu _{\mathrm {t} }{\bar {S}}_{ij}}" loading="lazy"></dd></dl>

where <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \nu _{\mathrm {t} }}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \nu _{\mathrm {t} }}</annotation>

</semantics>

</math><img src="./f3f2db1730054e31589bfae8d834d492ba613ad1.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:2.02ex; height:2.009ex;" alt="{\displaystyle \nu _{\mathrm {t} }}" loading="lazy"> is the turbulent eddy viscosity and <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\bar {S}}_{ij}={\frac {1}{2}}\left({\frac {\partial {\bar {u}}_{i}}{\partial x_{j}}}+{\frac {\partial {\bar {u}}_{j}}{\partial x_{i}}}\right)}">

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

<mfrac>

</mfrac>

</mrow>

<mrow>

<mrow>

<mfrac>

<mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mfrac>

<mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\bar {S}}_{ij}={\frac {1}{2}}\left({\frac {\partial {\bar {u}}_{i}}{\partial x_{j}}}+{\frac {\partial {\bar {u}}_{j}}{\partial x_{i}}}\right)}</annotation>

</semantics>

</math><img src="./0d2f2893c3d2a02601ebe177edd11d0fdbc8c34f.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.671ex; width:23.583ex; height:6.509ex;" alt="{\displaystyle {\bar {S}}_{ij}={\frac {1}{2}}\left({\frac {\partial {\bar {u}}_{i}}{\partial x_{j}}}+{\frac {\partial {\bar {u}}_{j}}{\partial x_{i}}}\right)}" loading="lazy"> is the rate-of-strain tensor.

Based on dimensional analysis, the eddy viscosity must have units of <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \left[\nu _{\mathrm {t} }\right]={\frac {\mathrm {m^{2}} }{\mathrm {s} }}}">

<mrow>

<msub>

</mrow>

</msub>

</mrow>

<mfrac>

<msup>

</mrow>

</msup>

</mrow>

</mrow>

</mfrac>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \left[\nu _{\mathrm {t} }\right]={\frac {\mathrm {m^{2}} }{\mathrm {s} }}}</annotation>

</semantics>

</math><img src="./184c1fccb884babf37ef4fde3f123842c969d309.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.838ex; width:10.239ex; height:5.676ex;" alt="{\displaystyle \left[\nu _{\mathrm {t} }\right]={\frac {\mathrm {m^{2}} }{\mathrm {s} }}}" loading="lazy">. Most eddy viscosity SGS models model the eddy viscosity as the product of a characteristic length scale and a characteristic velocity scale.

<div class="mw-heading mw-heading5"><h5 id="Smagorinsky–Lilly_model">Smagorinsky–Lilly model</h5></div>

The first SGS model developed was the Smagorinsky–Lilly SGS model, which was developed by <a href="Joseph_Smagorinsky" title="Joseph Smagorinsky">Smagorinsky</a><a href="#cite_note-Smagorinsky_1963-1">[1]</a> and used in the first LES simulation by Deardorff.<a href="#cite_note-Deardorff_1970-2">[2]</a> It models the eddy viscosity as:

<msub>

</mrow>

</msub>

<msup>

</mrow>

</msup>

<msqrt>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</msqrt>

</mrow>

<msup>

</mrow>

</msup>

<mrow>

<mover>

</mover>

</mrow>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \nu _{\mathrm {t} }=C\Delta ^{2}{\sqrt {2{\bar {S}}_{ij}{\bar {S}}_{ij}}}=C\Delta ^{2}\left|{\bar {S}}\right|}</annotation>

</semantics>

</math><img src="./37c3761fd42602de3ebacbf8284957b93e91c088.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.671ex; width:30.57ex; height:4.843ex;" alt="{\displaystyle \nu _{\mathrm {t} }=C\Delta ^{2}{\sqrt {2{\bar {S}}_{ij}{\bar {S}}_{ij}}}=C\Delta ^{2}\left|{\bar {S}}\right|}" loading="lazy"></dd></dl>

where <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \Delta }">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \Delta }</annotation>

</semantics>

</math><img src="./32769037c408874e1890f77554c65f39c523ebe2.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.936ex; height:2.176ex;" alt="{\displaystyle \Delta }" loading="lazy"> is the grid size and <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C}</annotation>

</semantics>

</math><img src="./4fc55753007cd3c18576f7933f6f089196732029.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.766ex; height:2.176ex;" alt="{\displaystyle C}" loading="lazy"> is a constant.

This method assumes that the energy production and dissipation of the small scales are in equilibrium - that is, <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \epsilon =\Pi }">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \epsilon =\Pi }</annotation>

</semantics>

</math><img src="./81b1a1035afa2425a09c764ecf6610353a2e72f0.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:5.786ex; height:2.176ex;" alt="{\displaystyle \epsilon =\Pi }" loading="lazy">.

<div class="mw-heading mw-heading5"><h5 id="The_Dynamic_Model_(Germano_et_al._and_beyond)">The Dynamic Model (Germano et al. and beyond)</h5></div>

Germano et al.<a href="#cite_note-Germano_1991-23">[23]</a> identified a number of studies using the Smagorinsky model that each found different values for the Smagorinsky constant <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C}</annotation>

</semantics>

</math><img src="./4fc55753007cd3c18576f7933f6f089196732029.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.766ex; height:2.176ex;" alt="{\displaystyle C}" loading="lazy"> for different flow configurations. In an attempt to formulate a more universal approach to SGS models, Germano et al. proposed a dynamic Smagorinsky model, which utilized two filters: a grid LES filter, denoted <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\overline {f}}}">

<mover>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\overline {f}}}</annotation>

</semantics>

</math><img src="./4de4554a5933ca51745edacd34a5a6aaea9a9f4e.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:1.484ex; height:3.343ex;" alt="{\displaystyle {\overline {f}}}" loading="lazy">, and a test LES filter, denoted <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\hat {f}}}">

<mover>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\hat {f}}}</annotation>

</semantics>

</math><img src="./14ce989fd75da938ec6f95a0cdb71037b23a11cb.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:1.699ex; height:3.176ex;" alt="{\displaystyle {\hat {f}}}" loading="lazy"> for any turbulent field <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle f}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle f}</annotation>

</semantics>

</math><img src="./132e57acb643253e7810ee9702d9581f159a1c61.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:1.279ex; height:2.509ex;" alt="{\displaystyle f}" loading="lazy">. The test filter is larger in size than the grid filter and adds an additional smoothing of the turbulence field over the already smoothed fields represented by the LES. Applying the test filter to the LES equations (which are obtained by applying the "grid" filter to Navier-Stokes equations) results in a new set of equations that are identical in form but with the SGS stress <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \tau _{ij}={\overline {u_{i}u_{j}}}-{\bar {u}}_{i}{\bar {u}}_{j}}">

<msub>

</mrow>

</msub>

<mover>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

</mover>

</mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \tau _{ij}={\overline {u_{i}u_{j}}}-{\bar {u}}_{i}{\bar {u}}_{j}}</annotation>

</semantics>

</math><img src="./15e2e14ac54b637d683e874dc63e74d15dbaef51.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:17.285ex; height:3.009ex;" alt="{\displaystyle \tau _{ij}={\overline {u_{i}u_{j}}}-{\bar {u}}_{i}{\bar {u}}_{j}}" loading="lazy"> replaced by <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle T_{ij}={\widehat {\overline {u_{i}u_{j}}}}-{\hat {\bar {u}}}_{i}{\hat {\bar {u}}}_{j}}">

<msub>

</mrow>

</msub>

<mover>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

</mover>

</mover>

</mrow>

<msub>

<mover>

<mover>

</mover>

</mrow>

</mover>

</mrow>

</mrow>

</msub>

<msub>

<mover>

<mover>

</mover>

</mrow>

</mover>

</mrow>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle T_{ij}={\widehat {\overline {u_{i}u_{j}}}}-{\hat {\bar {u}}}_{i}{\hat {\bar {u}}}_{j}}</annotation>

</semantics>

</math><img src="./e85159f25cdce2df68c1d02dade74a60981e4f2c.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:17.626ex; height:4.009ex;" alt="{\displaystyle T_{ij}={\widehat {\overline {u_{i}u_{j}}}}-{\hat {\bar {u}}}_{i}{\hat {\bar {u}}}_{j}}" loading="lazy">. Germano et al. noted that even though neither <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \tau _{ij}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \tau _{ij}}</annotation>

</semantics>

</math><img src="./97d0190e7517254a655e8033badbc01c8590381c.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:2.493ex; height:2.343ex;" alt="{\displaystyle \tau _{ij}}" loading="lazy"> nor <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle T_{ij}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle T_{ij}}</annotation>

</semantics>

</math><img src="./9de5760ded748359e36c7fb067c45f5e9642e890.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:2.835ex; height:2.843ex;" alt="{\displaystyle T_{ij}}" loading="lazy"> can be computed exactly because of the presence of unresolved scales, there is an exact relation connecting these two tensors. This relation, known as the Germano identity is

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle L_{ij}=T_{ij}-{\hat {\tau }}_{ij}.}</annotation>

</semantics>

</math><img src="./25c4718cca8cf9091b1a26c9f2fc5e52961c270d.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:15.325ex; height:2.843ex;" alt="{\displaystyle L_{ij}=T_{ij}-{\hat {\tau }}_{ij}.}" loading="lazy">

Here <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle L_{ij}={\widehat {{\bar {u}}_{i}{\bar {u}}_{j}}}-{\widehat {{\bar {u}}_{i}}}{\widehat {{\bar {u}}_{j}}}}">

<msub>

</mrow>

</msub>

<mover>

<mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mrow>

</mover>

</mrow>

<mover>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mover>

</mrow>

<mover>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle L_{ij}={\widehat {{\bar {u}}_{i}{\bar {u}}_{j}}}-{\widehat {{\bar {u}}_{i}}}{\widehat {{\bar {u}}_{j}}}}</annotation>

</semantics>

</math><img src="./838af6dbe7e27705843e22b8b1b84dd59268eb1c.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; margin-right: -0.009ex; width:18.044ex; height:3.509ex;" alt="{\displaystyle L_{ij}={\widehat {{\bar {u}}_{i}{\bar {u}}_{j}}}-{\widehat {{\bar {u}}_{i}}}{\widehat {{\bar {u}}_{j}}}}" loading="lazy"> can be explicitly evaluated as it involves only the filtered velocities and the operation of test filtering. The significance of the identity is that if one assumes that turbulence is self similar so that the SGS stress at the grid and test levels have the same form <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \tau _{ij}-(\tau _{kk}/3)\delta _{ij}=-2C\Delta ^{2}|{\bar {S}}_{ij}|{\bar {S}}_{ij}}">

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

<msub>

</mrow>

</msub>

<msup>

</mrow>

</msup>

</mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \tau _{ij}-(\tau _{kk}/3)\delta _{ij}=-2C\Delta ^{2}|{\bar {S}}_{ij}|{\bar {S}}_{ij}}</annotation>

</semantics>

</math><img src="./30b6a92e0d67fd9f0c0047013f6ae3bd10430629.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:33.158ex; height:3.343ex;" alt="{\displaystyle \tau _{ij}-(\tau _{kk}/3)\delta _{ij}=-2C\Delta ^{2}|{\bar {S}}_{ij}|{\bar {S}}_{ij}}" loading="lazy"> and <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle T_{ij}-(T_{kk}/3)\delta _{ij}=-2C{\hat {\Delta }}^{2}|{\hat {\bar {S}}}_{ij}|{\hat {\bar {S}}}_{ij}}">

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

<msub>

</mrow>

</msub>

<msup>

<mover>

</mover>

</mrow>

</mrow>

</msup>

</mrow>

<msub>

<mover>

<mover>

</mover>

</mrow>

</mover>

</mrow>

</mrow>

</msub>

</mrow>

<msub>

<mover>

<mover>

</mover>

</mrow>

</mover>

</mrow>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle T_{ij}-(T_{kk}/3)\delta _{ij}=-2C{\hat {\Delta }}^{2}|{\hat {\bar {S}}}_{ij}|{\hat {\bar {S}}}_{ij}}</annotation>

</semantics>

</math><img src="./9818b6c0c9da7666ca6b0c8bef93c6244e5b0d5f.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:33.84ex; height:4.009ex;" alt="{\displaystyle T_{ij}-(T_{kk}/3)\delta _{ij}=-2C{\hat {\Delta }}^{2}|{\hat {\bar {S}}}_{ij}|{\hat {\bar {S}}}_{ij}}" loading="lazy">, then the Germano identity provides an equation from which the Smagorinsky coefficient <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C}</annotation>

</semantics>

</math><img src="./4fc55753007cd3c18576f7933f6f089196732029.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.766ex; height:2.176ex;" alt="{\displaystyle C}" loading="lazy"> (which is no longer a 'constant') can potentially be determined.

[Inherent in the procedure is the assumption that the coefficient <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C}</annotation>

</semantics>

</math><img src="./4fc55753007cd3c18576f7933f6f089196732029.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.766ex; height:2.176ex;" alt="{\displaystyle C}" loading="lazy"> is invariant of scale (see review

<a href="#cite_note-MeneveauKatz_2000-24">[24]</a>)].

In order to do this, two additional steps were introduced in the original formulation. First, one assumed that even though <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C}</annotation>

</semantics>

</math><img src="./4fc55753007cd3c18576f7933f6f089196732029.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.766ex; height:2.176ex;" alt="{\displaystyle C}" loading="lazy"> was in principle variable, the variation was sufficiently slow that it can be moved out of the filtering operation <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\widehat {C(.)}}=C{\widehat {(.)}}}">

<mover>

<mrow>

</mrow>

</mover>

</mrow>

<mover>

<mrow>

</mrow>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\widehat {C(.)}}=C{\widehat {(.)}}}</annotation>

</semantics>

</math><img src="./69a572ebb88347028ef1c67ba70bf31b910030f5.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:12.318ex; height:3.676ex;" alt="{\displaystyle {\widehat {C(.)}}=C{\widehat {(.)}}}" loading="lazy">. Second, since <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C}</annotation>

</semantics>

</math><img src="./4fc55753007cd3c18576f7933f6f089196732029.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.766ex; height:2.176ex;" alt="{\displaystyle C}" loading="lazy"> was a scalar, the Germano identity was contracted with a second rank tensor (the rate of strain tensor was chosen) to convert it to a scalar equation from which <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C}</annotation>

</semantics>

</math><img src="./4fc55753007cd3c18576f7933f6f089196732029.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.766ex; height:2.176ex;" alt="{\displaystyle C}" loading="lazy"> could be determined.

Lilly

found a less arbitrary and therefore more satisfactory approach for obtaining C from the tensor identity. He noted that the Germano identity required the satisfaction of nine equations at each point in space (of which only five are independent) for a single quantity <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C}</annotation>

</semantics>

</math><img src="./4fc55753007cd3c18576f7933f6f089196732029.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.766ex; height:2.176ex;" alt="{\displaystyle C}" loading="lazy">. The problem of obtaining <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C}</annotation>

</semantics>

</math><img src="./4fc55753007cd3c18576f7933f6f089196732029.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.766ex; height:2.176ex;" alt="{\displaystyle C}" loading="lazy"> was therefore over-determined. He proposed therefore that <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C}</annotation>

</semantics>

</math><img src="./4fc55753007cd3c18576f7933f6f089196732029.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.766ex; height:2.176ex;" alt="{\displaystyle C}" loading="lazy"> be determined using a least square fit by minimizing the residuals. This results in

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C={\frac {L_{ij}m_{ij}}{m_{kl}m_{kl}}}.}</annotation>

</semantics>

</math><img src="./36a7131fef676982dc6d982929b3669faab35748.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.171ex; width:13.586ex; height:5.843ex;" alt="{\displaystyle C={\frac {L_{ij}m_{ij}}{m_{kl}m_{kl}}}.}" loading="lazy"></dd></dl>

Here

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle m_{ij}=\alpha _{ij}-{\widehat {\beta }}_{ij}}</annotation>

</semantics>

</math><img src="./8a3fc221d663ec74af2e8454564370e3105ff3c7.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.171ex; width:15.415ex; height:3.843ex;" alt="{\displaystyle m_{ij}=\alpha _{ij}-{\widehat {\beta }}_{ij}}" loading="lazy"></dd></dl>

and for brevity

<msub>

</mrow>

</msub>

<msup>

<mover>

</mover>

</mrow>

</mrow>

</msup>

</mrow>

<mover>

<mover>

</mover>

</mrow>

</mover>

</mrow>

</mrow>

<msub>

<mover>

<mover>

</mover>

</mrow>

</mover>

</mrow>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \alpha _{ij}=-2{\hat {\Delta }}^{2}|{\hat {\bar {S}}}|{\hat {\bar {S}}}_{ij}}</annotation>

</semantics>

</math><img src="./ddc59fd4939a856ffb0701023bcced0a464a8997.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:17.94ex; height:4.009ex;" alt="{\displaystyle \alpha _{ij}=-2{\hat {\Delta }}^{2}|{\hat {\bar {S}}}|{\hat {\bar {S}}}_{ij}}" loading="lazy">,

<msub>

</mrow>

</msub>

<msup>

</mrow>

</msup>

</mrow>

<mover>

</mover>

</mrow>

</mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \beta _{ij}=-2\Delta ^{2}|{\bar {S}}|{\bar {S}}_{ij}}</annotation>

</semantics>

</math><img src="./72576aab964f9cad38b7f74f3f68a926afef32df.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:17.768ex; height:3.343ex;" alt="{\displaystyle \beta _{ij}=-2\Delta ^{2}|{\bar {S}}|{\bar {S}}_{ij}}" loading="lazy">

Initial attempts to implement the model in LES simulations proved unsuccessful. First, the computed coefficient

was not at all "slowly varying" as assumed and varied as much as any other turbulent field. Secondly,

the computed <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C}</annotation>

</semantics>

</math><img src="./4fc55753007cd3c18576f7933f6f089196732029.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.766ex; height:2.176ex;" alt="{\displaystyle C}" loading="lazy"> could be positive as well as negative. The latter fact in itself should not be regarded as a

shortcoming as a priori tests using filtered DNS fields have shown that the local subgrid dissipation rate

<msub>

</mrow>

</msub>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle -\tau _{ij}{\bar {S}}_{ij}}</annotation>

</semantics>

</math><img src="./56176b9d8e44bda114077ca17adad8d155a7eebf.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:7.351ex; height:3.343ex;" alt="{\displaystyle -\tau _{ij}{\bar {S}}_{ij}}" loading="lazy"> in a turbulent field is almost as likely to be negative as it is positive even though the integral over the fluid domain is always positive representing a net dissipation of energy in the large scales. A slight preponderance of positive values as opposed to strict positivity of the eddy-viscosity results in the observed net dissipation. This so-called "backscatter" of energy from small to large scales indeed corresponds to negative C values in the Smagorinsky model. Nevertheless, the Germano-Lilly formulation was found not to result in stable calculations. An ad hoc measure was adopted by averaging the numerator and denominator over homogeneous directions (where such directions exist in the flow)

<mfrac>

<mrow>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

</mrow>

<mrow>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

</mrow>

</mfrac>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C={\frac {\left\langle L_{ij}m_{ij}\right\rangle }{\left\langle m_{kl}m_{kl}\right\rangle }}.}</annotation>

</semantics>

</math><img src="./cab67376125b551374337d7d96830715cc0113ee.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.671ex; width:15.396ex; height:6.509ex;" alt="{\displaystyle C={\frac {\left\langle L_{ij}m_{ij}\right\rangle }{\left\langle m_{kl}m_{kl}\right\rangle }}.}" loading="lazy"></dd></dl>

When the averaging involved a large enough statistical sample that the computed <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C}</annotation>

</semantics>

</math><img src="./4fc55753007cd3c18576f7933f6f089196732029.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.766ex; height:2.176ex;" alt="{\displaystyle C}" loading="lazy"> was positive (or at

least only rarely negative) stable calculations were possible. Simply setting the negative values to zero (a procedure called "clipping") with or without the averaging also resulted in stable calculations.

Meneveau proposed <a href="#cite_note-Meneveauetal_1996-26">[26]</a>

an averaging over Lagrangian fluid trajectories with an exponentially decaying "memory". This can be applied to problems lacking homogeneous directions and can be stable if the effective time over which the averaging is done is long enough and yet not so long as to smooth out spatial inhomogenieties of interest.

Lilly's modification of the Germano method followed by a statistical averaging or synthetic removal of negative viscosity regions seems ad hoc, even if it could be made to "work". An alternate formulation of the least square minimization procedure known as the "Dynamic Localization Model" (DLM) was suggested by

Ghosal et al.<a href="#cite_note-Ghosal_1995-27">[27]</a>

In this approach one first defines a quantity

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle E_{ij}=L_{ij}-T_{ij}+{\hat {\tau }}_{ij}}</annotation>

</semantics>

</math><img src="./96f18da78402504b21f713af76ffaf76fe15007e.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:20.711ex; height:2.843ex;" alt="{\displaystyle E_{ij}=L_{ij}-T_{ij}+{\hat {\tau }}_{ij}}" loading="lazy"></dd></dl>

with the tensors <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \tau _{ij}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \tau _{ij}}</annotation>

</semantics>

</math><img src="./97d0190e7517254a655e8033badbc01c8590381c.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:2.493ex; height:2.343ex;" alt="{\displaystyle \tau _{ij}}" loading="lazy"> and <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle T_{ij}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle T_{ij}}</annotation>

</semantics>

</math><img src="./9de5760ded748359e36c7fb067c45f5e9642e890.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:2.835ex; height:2.843ex;" alt="{\displaystyle T_{ij}}" loading="lazy"> replaced by the appropriate SGS model. This tensor then represents the amount by which the subgrid model fails to respect the Germano identity at each spatial location. In Lilly's approach, <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C}</annotation>

</semantics>

</math><img src="./4fc55753007cd3c18576f7933f6f089196732029.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.766ex; height:2.176ex;" alt="{\displaystyle C}" loading="lazy"> is then pulled out of the hat operator

<mover>

<mrow>

</mrow>

</mover>

</mrow>

<mover>

<mrow>

</mrow>

</mover>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\widehat {C(.)}}=C{\widehat {(.)}}}</annotation>

</semantics>

</math><img src="./69a572ebb88347028ef1c67ba70bf31b910030f5.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:12.318ex; height:3.676ex;" alt="{\displaystyle {\widehat {C(.)}}=C{\widehat {(.)}}}" loading="lazy"></dd></dl>

making <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle E_{ij}}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle E_{ij}}</annotation>

</semantics>

</math><img src="./acbcb625c128efafee881204113cd9e7a8a293a0.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:3.192ex; height:2.843ex;" alt="{\displaystyle E_{ij}}" loading="lazy"> an algebraic function of <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C}</annotation>

</semantics>

</math><img src="./4fc55753007cd3c18576f7933f6f089196732029.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.766ex; height:2.176ex;" alt="{\displaystyle C}" loading="lazy"> which is then determined by requiring that

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle E_{ij}E_{ij}}</annotation>

</semantics>

</math><img src="./27b3341c0bb4c8775cf13c898d77be712963d19a.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:6.385ex; height:2.843ex;" alt="{\displaystyle E_{ij}E_{ij}}" loading="lazy"> considered as a function of C have the least possible value.

However, since the <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C}</annotation>

</semantics>

</math><img src="./4fc55753007cd3c18576f7933f6f089196732029.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.766ex; height:2.176ex;" alt="{\displaystyle C}" loading="lazy"> thus obtained turns out to be just as variable as any other fluctuating quantity in turbulence, the original assumption of the constancy of <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C}</annotation>

</semantics>

</math><img src="./4fc55753007cd3c18576f7933f6f089196732029.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.766ex; height:2.176ex;" alt="{\displaystyle C}" loading="lazy"> cannot be justified a posteriori. In the DLM approach one avoids this inconsistency by not invoking the step of removing

C from the test filtering operation. Instead, one defines a global error over the entire flow domain by the quantity

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle E[C]=\int E_{ij}E_{ij}dV}</annotation>

</semantics>

</math><img src="./37685ef9c3d3b5c67410d4c450a7e4fa2526c695.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.338ex; width:19.903ex; height:5.676ex;" alt="{\displaystyle E[C]=\int E_{ij}E_{ij}dV}" loading="lazy"></dd></dl>

where the integral ranges over the whole fluid volume. This global error <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle E[C(x,y,z,t)]}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle E[C(x,y,z,t)]}</annotation>

</semantics>

</math><img src="./43664d34c18bf678de5159d42a59d37509ec5697.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:14.16ex; height:2.843ex;" alt="{\displaystyle E[C(x,y,z,t)]}" loading="lazy"> is then a functional of the spatially varying function <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C(x,y,z,t)}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C(x,y,z,t)}</annotation>

</semantics>

</math><img src="./4a105104ae37a639724b26c63b6faa00146b6b1f.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:11.09ex; height:2.843ex;" alt="{\displaystyle C(x,y,z,t)}" loading="lazy"> (here the time instant, <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle t}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle t}</annotation>

</semantics>

</math><img src="./65658b7b223af9e1acc877d848888ecdb4466560.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:0.84ex; height:2.009ex;" alt="{\displaystyle t}" loading="lazy">, is fixed and therefore appears just as a parameter) which is determined so as to minimize this functional. The solution to this variational problem is that <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C}</annotation>

</semantics>

</math><img src="./4fc55753007cd3c18576f7933f6f089196732029.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.766ex; height:2.176ex;" alt="{\displaystyle C}" loading="lazy"> must satisfy a Fredholm integral equation of the second kind

</mrow>

</mrow>

</mrow>

</mrow>

</mrow>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C({\boldsymbol {x}})=f({\boldsymbol {x}})+\int K({\boldsymbol {x}},{\boldsymbol {y}})C({\boldsymbol {y}})d{\boldsymbol {y}}}</annotation>

</semantics>

</math><img src="./87b5ecbe67039b19736065953a2eac9cba225168.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.338ex; width:33.594ex; height:5.676ex;" alt="{\displaystyle C({\boldsymbol {x}})=f({\boldsymbol {x}})+\int K({\boldsymbol {x}},{\boldsymbol {y}})C({\boldsymbol {y}})d{\boldsymbol {y}}}" loading="lazy"></dd></dl>

where the functions <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle K({\boldsymbol {x}},{\boldsymbol {y}})}">

</mrow>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle K({\boldsymbol {x}},{\boldsymbol {y}})}</annotation>

</semantics>

</math><img src="./4e6ab9b9ebe0924893607e6f0a0bcda5de6b273c.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:7.812ex; height:2.843ex;" alt="{\displaystyle K({\boldsymbol {x}},{\boldsymbol {y}})}" loading="lazy"> and <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle f({\boldsymbol {x}})}">

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle f({\boldsymbol {x}})}</annotation>

</semantics>

</math><img src="./63286363e334e9bedd72d688f8c3c3405aa1d939.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:4.62ex; height:2.843ex;" alt="{\displaystyle f({\boldsymbol {x}})}" loading="lazy"> are defined in terms of the resolved fields <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle L_{ij},\alpha _{ij},\beta _{ij}}">

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle L_{ij},\alpha _{ij},\beta _{ij}}</annotation>

</semantics>

</math><img src="./5d3207d6eb16d3a1dcda92c8dbbc7da8cafe9925.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:10.885ex; height:2.843ex;" alt="{\displaystyle L_{ij},\alpha _{ij},\beta _{ij}}" loading="lazy"> and are therefore known at each time step and the integral ranges over the whole fluid domain. The integral equation is solved numerically by an iteration procedure and convergence was found to be generally rapid if used with a pre-conditioning scheme. Even though this variational approach removes an inherent inconsistency in Lilly's approach, the <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C(x,y,z,t)}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C(x,y,z,t)}</annotation>

</semantics>

</math><img src="./4a105104ae37a639724b26c63b6faa00146b6b1f.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:11.09ex; height:2.843ex;" alt="{\displaystyle C(x,y,z,t)}" loading="lazy"> obtained from the integral equation still displayed the instability associated with negative viscosities. This can be resolved by insisting that <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle E[C]}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle E[C]}</annotation>

</semantics>

</math><img src="./70d52ebec4577f8617d3b2d70ccb501bf65f7797.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:4.836ex; height:2.843ex;" alt="{\displaystyle E[C]}" loading="lazy"> be minimized subject to the constraint <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C(x,y,z,t)\geq 0}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C(x,y,z,t)\geq 0}</annotation>

</semantics>

</math><img src="./2fc3260aa940cf14040b829cb2dc18bcbea0c16f.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:15.351ex; height:2.843ex;" alt="{\displaystyle C(x,y,z,t)\geq 0}" loading="lazy">. This leads to an equation for <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C}</annotation>

</semantics>

</math><img src="./4fc55753007cd3c18576f7933f6f089196732029.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.766ex; height:2.176ex;" alt="{\displaystyle C}" loading="lazy"> that is nonlinear

</mrow>

<msub>

<mrow>

<mrow>

</mrow>

</mrow>

</mrow>

</mrow>

</mrow>

</mrow>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C({\boldsymbol {x}})=\left[f({\boldsymbol {x}})+\int K({\boldsymbol {x}},{\boldsymbol {y}})C({\boldsymbol {y}})d{\boldsymbol {y}}\right]_{+}}</annotation>

</semantics>

</math><img src="./01fd8d503ee7dbcf432468ad08c669e7f31dd340.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.671ex; width:37.559ex; height:6.343ex;" alt="{\displaystyle C({\boldsymbol {x}})=\left[f({\boldsymbol {x}})+\int K({\boldsymbol {x}},{\boldsymbol {y}})C({\boldsymbol {y}})d{\boldsymbol {y}}\right]_{+}}" loading="lazy"></dd></dl>

Here the suffix + indicates the "positive part of" that is, <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle x_{+}=(x+|x|)/2}">

<msub>

</mrow>

</msub>

</mrow>

</mrow>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle x_{+}=(x+|x|)/2}</annotation>

</semantics>

</math><img src="./49ae90baa9a2ca5de36652e1511c4c5b5973039c.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:16.867ex; height:2.843ex;" alt="{\displaystyle x_{+}=(x+|x|)/2}" loading="lazy">. Even though this superficially looks like "clipping" it is not an ad hoc scheme but a bonafide solution of the constrained variational problem. This DLM(+) model was found to be stable and yielded excellent results for forced and decaying isotropic turbulence, channel flows and a variety of other more complex geometries. If a flow happens to have homogeneous directions (let us say the directions x and z) then one can introduce the ansatz

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C=C(y,t)}</annotation>

</semantics>

</math><img src="./59defd6aae18a169948294cd25fbebc018e2ca31.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:11.469ex; height:2.843ex;" alt="{\displaystyle C=C(y,t)}" loading="lazy">. The variational approach then immediately yields Lilly's result with averaging over homogeneous directions without any need for ad hoc modifications of a prior result.

One shortcoming of the DLM(+) model was that it did not describe backscatter which is known to be a real "thing" from analyzing DNS data. Two approaches were developed to address this. In one approach due to Carati et al.

a fluctuating force with amplitude determined by the fluctuation-dissipation theorem is added in

analogy to Landau's theory of fluctuating hydrodynamics. In the second approach, one notes that

any "backscattered" energy appears in the resolved scales only at the expense of energy in the subgrid

scales. The DLM can be modified in a simple way to take into account this physical fact so as to allow

for backscatter while being inherently stable. This k-equation version of the DLM, DLM(k) replaces

</mrow>

<mover>

</mover>

</mrow>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \Delta |{\bar {S}}|}</annotation>

</semantics>

</math><img src="./1e45c2d390cda8015e9cf118d805e49096e9ccfa.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:4.802ex; height:3.176ex;" alt="{\displaystyle \Delta |{\bar {S}}|}" loading="lazy"> in the Smagorinsky eddy viscosity model by <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle {\sqrt {k}}}">

<msqrt>

</msqrt>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\sqrt {k}}}</annotation>

</semantics>

</math><img src="./de05da7afd02cfd22d059bf17a0c165e3079d5d4.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:3.147ex; height:3.009ex;" alt="{\displaystyle {\sqrt {k}}}" loading="lazy"> as an appropriate velocity scale. The procedure for determining <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C}">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C}</annotation>

</semantics>

</math><img src="./4fc55753007cd3c18576f7933f6f089196732029.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.766ex; height:2.176ex;" alt="{\displaystyle C}" loading="lazy"> remains identical to the "unconstrained" version except that the tensors <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \alpha _{ij}=-2{\hat {\Delta }}{\sqrt {K}}{\hat {\bar {S}}}_{ij}}">

<msub>

</mrow>

</msub>

<mover>

</mover>

</mrow>

<msqrt>

</msqrt>

</mrow>

<msub>

<mover>

<mover>

</mover>

</mrow>

</mover>

</mrow>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \alpha _{ij}=-2{\hat {\Delta }}{\sqrt {K}}{\hat {\bar {S}}}_{ij}}</annotation>

</semantics>

</math><img src="./ed24552d87e0c04bce5aa5e679c5cdd03f50a87d.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:18.021ex; height:3.843ex;" alt="{\displaystyle \alpha _{ij}=-2{\hat {\Delta }}{\sqrt {K}}{\hat {\bar {S}}}_{ij}}" loading="lazy">,

<msub>

</mrow>

</msub>

<mover>

</mover>

</mrow>

<msqrt>

</msqrt>

</mrow>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \beta _{ij}=-2{\hat {\Delta }}{\sqrt {k}}{\bar {S}}_{ij}}</annotation>

</semantics>

</math><img src="./eae4e2c907c5e909cc32b6d05a9563742447214e.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -1.005ex; width:16.995ex; height:3.509ex;" alt="{\displaystyle \beta _{ij}=-2{\hat {\Delta }}{\sqrt {k}}{\bar {S}}_{ij}}" loading="lazy"> where the sub-test scale kinetic

energy K is related to the subgrid scale kinetic energy k by <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle K=k+L_{ii}/2}">

<msub>

</mrow>

</msub>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle K=k+L_{ii}/2}</annotation>

</semantics>

</math><img src="./d8522741a9507c4d604dfafdb15a4d03a544082b.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.838ex; width:14.491ex; height:2.843ex;" alt="{\displaystyle K=k+L_{ii}/2}" loading="lazy">

(follows by taking the trace of the Germano identity). To determine k we now use a transport equation

<mfrac>

<mrow>

</mrow>

<mrow>

</mrow>

</mfrac>

</mrow>

<msub>

</mrow>

</msub>

<mfrac>

<mrow>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<msub>

</mrow>

</msub>

<msub>

<mover>

</mover>

</mrow>

</mrow>

</msub>

<mfrac>

<msub>

</mrow>

</msub>

</mfrac>

</mrow>

<msup>

</mrow>

</mrow>

</msup>

<mfrac>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

<mrow>

<mrow>

<msqrt>

</msqrt>

</mrow>

<mfrac>

<mrow>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

</mrow>

<mfrac>

<mrow>

<msup>

</mrow>

</msup>

</mrow>

<mrow>

<msub>

</mrow>

</msub>

<msub>

</mrow>

</msub>

</mrow>

</mfrac>

</mrow>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle {\frac {\partial k}{\partial t}}+u_{j}{\frac {\partial k}{\partial x_{j}}}=-\tau _{ij}{\bar {S}}_{ij}-{\frac {C_{*}}{\Delta }}k^{3/2}+{\frac {\partial }{\partial x_{j}}}\left(D\Delta {\sqrt {k}}{\frac {\partial k}{\partial x_{j}}}\right)+\nu {\frac {\partial ^{2}k}{\partial x_{j}\partial x_{j}}}}</annotation>

</semantics>

</math><img src="./76d43f13149070a2625c07453345907511e875e9.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -2.671ex; width:68.057ex; height:6.509ex;" alt="{\displaystyle {\frac {\partial k}{\partial t}}+u_{j}{\frac {\partial k}{\partial x_{j}}}=-\tau _{ij}{\bar {S}}_{ij}-{\frac {C_{*}}{\Delta }}k^{3/2}+{\frac {\partial }{\partial x_{j}}}\left(D\Delta {\sqrt {k}}{\frac {\partial k}{\partial x_{j}}}\right)+\nu {\frac {\partial ^{2}k}{\partial x_{j}\partial x_{j}}}}" loading="lazy"></dd></dl>

where <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle \nu }">

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle \nu }</annotation>

</semantics>

</math><img src="./c15bbbb971240cf328aba572178f091684585468.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.338ex; width:1.232ex; height:1.676ex;" alt="{\displaystyle \nu }" loading="lazy"> is the kinematic viscosity and <math xmlns="http://www.w3.org/1998/Math/MathML" alttext="{\displaystyle C_{*},D}">

<msub>

</mrow>

</msub>

</mstyle>

</mrow>

<annotation encoding="application/x-tex">{\displaystyle C_{*},D}</annotation>

</semantics>

</math><img src="./95f65fb5cd1ac092e29193c092215a6f99cb41f0.svg" class="mwe-math-fallback-image-inline mw-invert skin-invert" aria-hidden="true" style="vertical-align: -0.671ex; width:5.674ex; height:2.509ex;" alt="{\displaystyle C_{*},D}" loading="lazy"> are positive coefficients

representing kinetic energy dissipation and diffusion respectively. These can be determined following the dynamic

procedure with constrained minimization as in DLM(+). This approach, though more expensive to implement than the DLM(+) was found to be stable and resulted in good agreement with experimental data for a variety of flows

tested. Furthermore, it is mathematically impossible for the DLM(k) to result in an unstable computation as the sum of the large scale and SGS energies is non-increasing by construction. Both of these approaches incorporating backscatter works well. They yield models that are slightly less dissipative with somewhat improved performance over the DLM(+). The DLM(k) model additionally yields the subgrid kinetic energy, which may be a physical quantity of interest. These improvements are achieved at a somewhat increased cost in model implementation.

The Dynamic Model originated at the 1990 <a rel="nofollow" class="external text" href="https://ctr.stanford.edu/ctr-summer-program">Summer Program</a> of the <a href="Center_for_Turbulence_Research" title="Center for Turbulence Research">Center for Turbulence Research</a> (CTR) at <a href="Stanford_University" title="Stanford University">Stanford University</a>. A series of "CTR-Tea" seminars celebrated the <a rel="nofollow" class="external text" href="https://ctr.stanford.edu/event-type/ctr-tea">30th Anniversary</a> <a rel="nofollow" class="external text" href="https://web.archive.org/web/20221030050059/https://ctr.stanford.edu/event-type/ctr-tea">Archived</a> 2022-10-30 at the <a href="Wayback_Machine" title="Wayback Machine">Wayback Machine</a> of this important milestone in turbulence modeling.

<div class="mw-heading mw-heading4"><h4 id="Structural_models">Structural models</h4></div>

/* start https://en.wikipedia.org/ */

.mw-parser-output .ambox{border:1px solid #a2a9b1;border-left:10px solid #36c;background-color:#fbfbfb;box-sizing:border-box}.mw-parser-output .ambox+link+.ambox,.mw-parser-output .ambox+link+style+.ambox,.mw-parser-output .ambox+link+link+.ambox,.mw-parser-output .ambox+.mw-empty-elt+link+.ambox,.mw-parser-output .ambox+.mw-empty-elt+link+style+.ambox,.mw-parser-output .ambox+.mw-empty-elt+link+link+.ambox{margin-top:-1px}html body.mediawiki .mw-parser-output .ambox.mbox-small-left{margin:4px 1em 4px 0;overflow:hidden;width:238px;border-collapse:collapse;font-size:88%;line-height:1.25em}.mw-parser-output .ambox-speedy{border-left:10px solid #b32424;background-color:#fee7e6}.mw-parser-output .ambox-delete{border-left:10px solid #b32424}.mw-parser-output .ambox-content{border-left:10px solid #f28500}.mw-parser-output .ambox-style{border-left:10px solid #fc3}.mw-parser-output .ambox-move{border-left:10px solid #9932cc}.mw-parser-output .ambox-protection{border-left:10px solid #a2a9b1}.mw-parser-output .ambox .mbox-text{border:none;padding:0.25em 0.5em;width:100%}.mw-parser-output .ambox .mbox-image{border:none;padding:2px 0 2px 0.5em;text-align:center}.mw-parser-output .ambox .mbox-imageright{border:none;padding:2px 0.5em 2px 0;text-align:center}.mw-parser-output .ambox .mbox-empty-cell{border:none;padding:0;width:1px}.mw-parser-output .ambox .mbox-image-div{width:52px}@media(min-width:720px){.mw-parser-output .ambox{margin:0 10%}}@media print{body.ns-0 .mw-parser-output .ambox{display:none!important}}

/* end https://en.wikipedia.org/ */

</style>

<ul><li><a href="Direct_numerical_simulation" title="Direct numerical simulation">Direct numerical simulation</a></li>

<li><a href="Fluid_mechanics" title="Fluid mechanics">Fluid mechanics</a></li>

<li><a href="Galilean_invariance" title="Galilean invariance">Galilean invariance</a> – an important property of certain types of filters</li>

<li><a href="Reynolds-averaged_Navier%E2%80%93Stokes_equations" title="Reynolds-averaged Navier–Stokes equations">Reynolds-averaged Navier–Stokes equations</a></li>

<li><a href="Turbulence" title="Turbulence">Turbulence</a></li></ul>

<div class="mw-heading mw-heading2"><h2 id="Further_reading">Further reading</h2></div>

<ul><li>Heus, T.; van Heerwaarden, C. C.;Jonker, H. J. J.; Pier Siebesma, A.; Axelsen, S.«<a rel="nofollow" class="external text" href="https://www.geosci-model-dev.net/3/415/2010/gmd-3-415-2010.html">Formulation of the Dutch Atmospheric Large-Eddy Simulation (DALES) and overview of its applications</a>» Geoscientific Model Development, 3, 2, 30-09-2010, pàg. 415–444. <a href="https://ca.wikipedia.org/wiki/DOI" class="extiw external" title="ca:DOI">DOI</a>: <a rel="nofollow" class="external text" href="https://www.geosci-model-dev.net/3/415/2010/gmd-3-415-2010.html">10.5194/gmd-3-415-2010</a>. <a href="https://ca.wikipedia.org/wiki/ISSN" class="extiw external" title="ca:ISSN">ISSN</a>: 1991-9603.</li></ul>

<div class="mw-heading mw-heading2"><h2 id="References">References</h2></div>

/* start https://en.wikipedia.org/ */

.mw-parser-output .reflist{margin-bottom:0.5em;list-style-type:decimal}@media screen{.mw-parser-output .reflist{font-size:90%}}.mw-parser-output .reflist .references{font-size:100%;margin-bottom:0;list-style-type:inherit}.mw-parser-output .reflist-columns-2{column-width:30em}.mw-parser-output .reflist-columns-3{column-width:25em}.mw-parser-output .reflist-columns{margin-top:0.3em}.mw-parser-output .reflist-columns ol{margin-top:0}.mw-parser-output .reflist-columns li{page-break-inside:avoid;break-inside:avoid-column}.mw-parser-output .reflist-upper-alpha{list-style-type:upper-alpha}.mw-parser-output .reflist-upper-roman{list-style-type:upper-roman}.mw-parser-output .reflist-lower-alpha{list-style-type:lower-alpha}.mw-parser-output .reflist-lower-greek{list-style-type:lower-greek}.mw-parser-output .reflist-lower-roman{list-style-type:lower-roman}

/* end https://en.wikipedia.org/ */

</style><div class="reflist">

/* start https://en.wikipedia.org/ */

.mw-parser-output cite.citation{font-style:inherit;word-wrap:break-word}.mw-parser-output .citation q{quotes:"\"""\"""'""'"}.mw-parser-output .citation:target{background-color:rgba(0,127,255,0.133)}.mw-parser-output .id-lock-free.id-lock-free a{background:url("./mw/Lock-green.svg")right 0.1em center/9px no-repeat}.mw-parser-output .id-lock-limited.id-lock-limited a,.mw-parser-output .id-lock-registration.id-lock-registration a{background:url("./mw/Lock-gray-alt-2.svg")right 0.1em center/9px no-repeat}.mw-parser-output .id-lock-subscription.id-lock-subscription a{background:url("./mw/Lock-red-alt-2.svg")right 0.1em center/9px no-repeat}.mw-parser-output .cs1-ws-icon a{background:url("./mw/Wikisource-logo.svg")right 0.1em center/12px no-repeat}body:not(.skin-timeless):not(.skin-minerva) .mw-parser-output .id-lock-free a,body:not(.skin-timeless):not(.skin-minerva) .mw-parser-output .id-lock-limited a,body:not(.skin-timeless):not(.skin-minerva) .mw-parser-output .id-lock-registration a,body:not(.skin-timeless):not(.skin-minerva) .mw-parser-output .id-lock-subscription a,body:not(.skin-timeless):not(.skin-minerva) .mw-parser-output .cs1-ws-icon a{background-size:contain;padding:0 1em 0 0}.mw-parser-output .cs1-code{color:inherit;background:inherit;border:none;padding:inherit}.mw-parser-output .cs1-hidden-error{display:none;color:var(--color-error,#d33)}.mw-parser-output .cs1-visible-error{color:var(--color-error,#d33)}.mw-parser-output .cs1-maint{display:none;color:#085;margin-left:0.3em}.mw-parser-output .cs1-kern-left{padding-left:0.2em}.mw-parser-output .cs1-kern-right{padding-right:0.2em}.mw-parser-output .citation .mw-selflink{font-weight:inherit}@media screen{.mw-parser-output .cs1-format{font-size:95%}html.skin-theme-clientpref-night .mw-parser-output .cs1-maint{color:#18911f}}@media screen and (prefers-color-scheme:dark){html.skin-theme-clientpref-os .mw-parser-output .cs1-maint{color:#18911f}}

/* end https://en.wikipedia.org/ */

</style><cite id="CITEREFSmagorinsky1963" class="citation journal cs1">Smagorinsky, Joseph (March 1963). <a rel="nofollow" class="external text" href="https://doi.org/10.1175%2F1520-0493%281963%29091%3C0099%3AGCEWTP%3E2.3.CO%3B2">"General Circulation Experiments with the Primitive Equations"</a>. Monthly Weather Review. 91 (3): 99–164. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/1963MWRv...91...99S">1963MWRv...91...99S</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1175%2F1520-0493%281963%29091%3C0099%3AGCEWTP%3E2.3.CO%3B2">10.1175/1520-0493(1963)091<0099:GCEWTP>2.3.CO;2</a>.</cite>

</li>

<li id="cite_note-Deardorff_1970-2">^ <a href="#cite_ref-Deardorff_1970_2-0">a</a> <a href="#cite_ref-Deardorff_1970_2-1">b</a> <cite id="CITEREFDeardorff1970" class="citation journal cs1">Deardorff, James (1970). "A numerical study of three-dimensional turbulent channel flow at large Reynolds numbers". <a href="Journal_of_Fluid_Mechanics" title="Journal of Fluid Mechanics">Journal of Fluid Mechanics</a>. 41 (2): 453–480. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/1970JFM....41..453D">1970JFM....41..453D</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1017%2FS0022112070000691">10.1017/S0022112070000691</a>. <a href="S2CID_(identifier)" class="mw-redirect" title="S2CID (identifier)">S2CID</a> <a rel="nofollow" class="external text" href="https://api.semanticscholar.org/CorpusID:121884175">121884175</a>.</cite>

</li>

<li id="cite_note-Pitsch_2006-3">^ <a href="#cite_ref-Pitsch_2006_3-0">a</a> <a href="#cite_ref-Pitsch_2006_3-1">b</a> <cite id="CITEREFPitsch2006" class="citation journal cs1">Pitsch, Heinz (2006). <a rel="nofollow" class="external text" href="http://caltechconf.library.caltech.edu/145/1/Abstract1.pdf">"Large-Eddy Simulation of Turbulent Combustion"</a> (PDF). Annual Review of Fluid Mechanics. 38 (1): 453–482. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/2006AnRFM..38..453P">2006AnRFM..38..453P</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1146%2Fannurev.fluid.38.050304.092133">10.1146/annurev.fluid.38.050304.092133</a>. <a href="S2CID_(identifier)" class="mw-redirect" title="S2CID (identifier)">S2CID</a> <a rel="nofollow" class="external text" href="https://api.semanticscholar.org/CorpusID:5487815">5487815</a>.</cite>

</li>

<li id="cite_note-Wagner_2007-4"><a href="#cite_ref-Wagner_2007_4-0">^</a> <cite id="CITEREFWagnerHüttlSagaut2007" class="citation book cs1">Wagner, Claus; Hüttl, Thomas; Sagaut, Pierre (2007). Large-Eddy Simulation for Acoustics. Cambridge University Press. <a href="ISBN_(identifier)" class="mw-redirect" title="ISBN (identifier)">ISBN</a> <bdi>978-0-521-87144-0</bdi>.</cite>

</li>

<li id="cite_note-Sullivan_1994-5"><a href="#cite_ref-Sullivan_1994_5-0">^</a> <cite id="CITEREFSullivanMcWilliamsMoeng1994" class="citation journal cs1">Sullivan, Peter P.; McWilliams, James C.; Moeng, Chin-Hoh (1994). "A subgrid-scale model for large-eddy simulation of planetary boundary-layer flows". Boundary-Layer Meteorology. 71 (3): 247–276. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/1994BoLMe..71..247S">1994BoLMe..71..247S</a>. <a href="CiteSeerX_(identifier)" class="mw-redirect" title="CiteSeerX (identifier)">CiteSeerX</a> <a rel="nofollow" class="external text" href="https://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.463.6006">10.1.1.463.6006</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1007%2FBF00713741">10.1007/BF00713741</a>. <a href="ISSN_(identifier)" class="mw-redirect" title="ISSN (identifier)">ISSN</a> <a rel="nofollow" class="external text" href="https://search.worldcat.org/issn/0006-8314">0006-8314</a>. <a href="S2CID_(identifier)" class="mw-redirect" title="S2CID (identifier)">S2CID</a> <a rel="nofollow" class="external text" href="https://api.semanticscholar.org/CorpusID:53051046">53051046</a>.</cite>

</li>

<li id="cite_note-piomelli-02-6"><a href="#cite_ref-piomelli-02_6-0">^</a> <cite id="CITEREFPiomelli,_UgoElias_Balaras2002" class="citation journal cs1">Piomelli, Ugo; Elias Balaras (2002). "Wall-layer models for large-eddy simulations". Annual Review of Fluid Mechanics. 34 (34): 349–374. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/2002AnRFM..34..349P">2002AnRFM..34..349P</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1146%2Fannurev.fluid.34.082901.144919">10.1146/annurev.fluid.34.082901.144919</a>.</cite>

</li>

<li id="cite_note-spalart-09-7"><a href="#cite_ref-spalart-09_7-0">^</a> <cite id="CITEREFSpalart2009" class="citation journal cs1">Spalart, Philippe R. (2009). "Detached-eddy simulation". Annual Review of Fluid Mechanics. 41 (1): 181–202. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/2009AnRFM..41..181S">2009AnRFM..41..181S</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1146%2Fannurev.fluid.010908.165130">10.1146/annurev.fluid.010908.165130</a>.</cite>

</li>

<li id="cite_note-fox-12-8"><a href="#cite_ref-fox-12_8-0">^</a> <cite id="CITEREFFox,_R._O.2012" class="citation journal cs1">Fox, R. O. (2012). "Large-eddy-simulation tools for multiphase flows". Annual Review of Fluid Mechanics. 44 (1): 47–76. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/2012AnRFM..44...47F">2012AnRFM..44...47F</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1146%2Fannurev-fluid-120710-101118">10.1146/annurev-fluid-120710-101118</a>.</cite>

</li>

<li id="cite_note-Pope_2000-9">^ <a href="#cite_ref-Pope_2000_9-0">a</a> <a href="#cite_ref-Pope_2000_9-1">b</a> <a href="#cite_ref-Pope_2000_9-2">c</a> <a href="#cite_ref-Pope_2000_9-3">d</a> <cite id="CITEREFPope,_S._B.2000" class="citation book cs1">Pope, S. B. (2000). Turbulent Flows. Cambridge University Press.</cite>

</li>

<li id="cite_note-Sagaut_2006-10">^ <a href="#cite_ref-Sagaut_2006_10-0">a</a> <a href="#cite_ref-Sagaut_2006_10-1">b</a> <cite id="CITEREFSagaut,_Pierre2006" class="citation book cs1">Sagaut, Pierre (2006). Large Eddy Simulation for Incompressible Flows (Third ed.). Springer. <a href="ISBN_(identifier)" class="mw-redirect" title="ISBN (identifier)">ISBN</a> <bdi>978-3-540-26344-9</bdi>.</cite>

</li>

<li id="cite_note-Leonard_1974-11">^ <a href="#cite_ref-Leonard_1974_11-0">a</a> <a href="#cite_ref-Leonard_1974_11-1">b</a> <cite id="CITEREFLeonard1975" class="citation book cs1">Leonard, A. (1975). "Energy cascade in large-eddy simulations of turbulent fluid flows". Turbulent Diffusion in Environmental Pollution, Proceedings of a Symposium held at Charlottesville. Advances in Geophysics A. Vol. 18. pp. 237–248. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/1975AdGeo..18..237L">1975AdGeo..18..237L</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1016%2FS0065-2687%2808%2960464-1">10.1016/S0065-2687(08)60464-1</a>. <a href="ISBN_(identifier)" class="mw-redirect" title="ISBN (identifier)">ISBN</a> <bdi>9780120188185</bdi>.</cite>

</li>

<li id="cite_note-Clark-12"><a href="#cite_ref-Clark_12-0">^</a> <cite id="CITEREFClarkFerzigerReynolds1979" class="citation journal cs1">Clark, R.; Ferziger, J.; <a href="William_Craig_Reynolds" title="William Craig Reynolds">Reynolds, W.</a> (1979). "Evaluation of subgrid-scale models using an accurately simulated turbulent flow". <a href="Journal_of_Fluid_Mechanics" title="Journal of Fluid Mechanics">Journal of Fluid Mechanics</a>. 91: 1–16. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/1979JFM....91....1C">1979JFM....91....1C</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1017%2FS002211207900001X">10.1017/S002211207900001X</a>. <a href="S2CID_(identifier)" class="mw-redirect" title="S2CID (identifier)">S2CID</a> <a rel="nofollow" class="external text" href="https://api.semanticscholar.org/CorpusID:120228458">120228458</a>.</cite>

</li>

<li id="cite_note-Favre_1983-13"><a href="#cite_ref-Favre_1983_13-0">^</a> <cite id="CITEREFFavre,_Alexandre1983" class="citation journal cs1">Favre, Alexandre (1983). "Turbulence: space-time statistical properties and behavior in supersonic flows". Physics of Fluids A. 23 (10): 2851–2863. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/1983PhFl...26.2851F">1983PhFl...26.2851F</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1063%2F1.864049">10.1063/1.864049</a>.</cite>

</li>

<li id="cite_note-Vreman_1995-14"><a href="#cite_ref-Vreman_1995_14-0">^</a> <cite id="CITEREFVreman,_BertGeurts,_BernardKuerten,_Hans1995" class="citation journal cs1">Vreman, Bert; Geurts, Bernard; Kuerten, Hans (1995). <a rel="nofollow" class="external text" href="https://research.utwente.nl/en/publications/subgridmodelling-in-les-of-compressible-flow(7c54958d-ebdd-4422-bf9d-c17052984a68).html">"Subgrid-modelling in LES of compressible flow"</a>. <a href="Applied_Scientific_Research" class="mw-redirect" title="Applied Scientific Research">Applied Scientific Research</a>. 45 (3): 191–203. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/1995FTC....54..191V">1995FTC....54..191V</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1007%2FBF00849116">10.1007/BF00849116</a>.</cite>

</li>

<li id="cite_note-Sagaut_2009-15"><a href="#cite_ref-Sagaut_2009_15-0">^</a> <cite id="CITEREFGarnier,_E.Adams,_N.Sagaut,_P.2009" class="citation book cs1">Garnier, E.; Adams, N.; Sagaut, P. (2009). <a rel="nofollow" class="external text" href="https://cds.cern.ch/record/1339029">Large eddy simulation for compressible flows</a>. Springer. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1007%2F978-90-481-2819-8">10.1007/978-90-481-2819-8</a>. <a href="ISBN_(identifier)" class="mw-redirect" title="ISBN (identifier)">ISBN</a> <bdi>978-90-481-2818-1</bdi>.</cite>

</li>

<li id="cite_note-Piomelli_1991-16"><a href="#cite_ref-Piomelli_1991_16-0">^</a> <cite id="CITEREFPiomelli,_U.Cabot,_W.Moin,_P.Lee,_S.1991" class="citation journal cs1">Piomelli, U.; Cabot, W.; <a href="Parviz_Moin" title="Parviz Moin">Moin, P.</a>; Lee, S. (1991). "Subgrid-scale backscatter in turbulent and transitional flows". Physics of Fluids A. 3 (7): 1766–1771. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/1991PhFlA...3.1766P">1991PhFlA...3.1766P</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1063%2F1.857956">10.1063/1.857956</a>. <a href="S2CID_(identifier)" class="mw-redirect" title="S2CID (identifier)">S2CID</a> <a rel="nofollow" class="external text" href="https://api.semanticscholar.org/CorpusID:54904570">54904570</a>.</cite>

</li>

<li id="cite_note-Ghosal_1996-17"><a href="#cite_ref-Ghosal_1996_17-0">^</a> <cite id="CITEREFGhosal,_S.1996" class="citation journal cs1">Ghosal, S. (April 1996). <a rel="nofollow" class="external text" href="https://doi.org/10.1006%2Fjcph.1996.0088">"An analysis of numerical errors in large-eddy simulations of turbulence"</a>. <a href="Journal_of_Computational_Physics" title="Journal of Computational Physics">Journal of Computational Physics</a>. 125 (1): 187–206. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/1996JCoPh.125..187G">1996JCoPh.125..187G</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1006%2Fjcph.1996.0088">10.1006/jcph.1996.0088</a>.</cite>

</li>

<li id="cite_note-Leveque_1992-18"><a href="#cite_ref-Leveque_1992_18-0">^</a> <cite id="CITEREFRandall_J._Leveque1992" class="citation book cs1">Randall J. Leveque (1992). Numerical Methods for Conservation Laws (2nd ed.). Birkhäuser Basel. <a href="ISBN_(identifier)" class="mw-redirect" title="ISBN (identifier)">ISBN</a> <bdi>978-3-7643-2723-1</bdi>.</cite>

</li>

<li id="cite_note-Grinstein_2007-19"><a href="#cite_ref-Grinstein_2007_19-0">^</a> <cite id="CITEREFGrinstein,_FernandoMargolin,_LenRider,_William2007" class="citation book cs1">Grinstein, Fernando; Margolin, Len; Rider, William (2007). Implicit large eddy simulation. Cambridge University Press. <a href="ISBN_(identifier)" class="mw-redirect" title="ISBN (identifier)">ISBN</a> <bdi>978-0-521-86982-9</bdi>.</cite>

</li>

<li id="cite_note-Li_P-20">^ <a href="#cite_ref-Li_P_20-0">a</a> <a href="#cite_ref-Li_P_20-1">b</a> <a href="#cite_ref-Li_P_20-2">c</a> <cite id="CITEREFLiEckelsMannZhang2018" class="citation journal cs1">Li, P.; Eckels, S.; Mann, G.; Zhang, N. (2018). "A Method of Measuring Turbulent Flow Structures With Particle Image Velocimetry and Incorporating Into Boundary Conditions of Large Eddy Simulations". Journal of Fluids Engineering. 140 (7). ASME International. 071401-071401-11. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1115%2F1.4039256">10.1115/1.4039256</a>.</cite>

</li>

<li id="cite_note-21"><a href="#cite_ref-21">^</a> <cite id="CITEREFTaborBaba-Ahmadi2010" class="citation journal cs1">Tabor, G. R.; Baba-Ahmadi, M. H. (2010). "Inlet conditions for large eddy simulation: a review". Computers & Fluids. 39 (4): 553–567. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1016%2Fj.compfluid.2009.10.007">10.1016/j.compfluid.2009.10.007</a>.</cite>

</li>

<li id="cite_note-Meneveau_2010-22"><a href="#cite_ref-Meneveau_2010_22-0">^</a> <cite id="CITEREFMeneveau,_C.2010" class="citation journal cs1">Meneveau, C. (2010). <a rel="nofollow" class="external text" href="https://doi.org/10.4249%2Fscholarpedia.9489">"Turbulence: Subgrid-Scale Modeling"</a>. Scholarpedia. 5 (1): 9489. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/2010SchpJ...5.9489M">2010SchpJ...5.9489M</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.4249%2Fscholarpedia.9489">10.4249/scholarpedia.9489</a>.</cite>

</li>

<li id="cite_note-Germano_1991-23"><a href="#cite_ref-Germano_1991_23-0">^</a> <cite id="CITEREFGermano,_M.Piomelli,_U.Moin,_P.Cabot,_W.1991" class="citation journal cs1">Germano, M.; Piomelli, U.; <a href="Parviz_Moin" title="Parviz Moin">Moin, P.</a>; Cabot, W. (1991). "A dynamic subgrid-scale eddy viscosity model". <a href="Physics_of_Fluids" title="Physics of Fluids">Physics of Fluids A</a>. 3 (7): 1760–1765. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/1991PhFlA...3.1760G">1991PhFlA...3.1760G</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1063%2F1.857955">10.1063/1.857955</a>. <a href="S2CID_(identifier)" class="mw-redirect" title="S2CID (identifier)">S2CID</a> <a rel="nofollow" class="external text" href="https://api.semanticscholar.org/CorpusID:55719851">55719851</a>.</cite>

</li>

<li id="cite_note-MeneveauKatz_2000-24"><a href="#cite_ref-MeneveauKatz_2000_24-0">^</a> <cite id="CITEREFMeneveau,_C.Katz,_J.2000" class="citation journal cs1">Meneveau, C.; Katz, J. (2000). "Scale-Invariance and Turbulence Models for Large-Eddy Simulation". Annu. Rev. Fluid Mech. 32 (1): 1–32. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/2000AnRFM..32....1M">2000AnRFM..32....1M</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1146%2Fannurev.fluid.32.1.1">10.1146/annurev.fluid.32.1.1</a>.</cite>

</li>

<li id="cite_note-Lilly_1992-25"><a href="#cite_ref-Lilly_1992_25-0">^</a> <cite id="CITEREFLilly1992" class="citation journal cs1">Lilly, D. K. (1992). "A proposed modification of the Germano subgrid-scale closure method". Physics of Fluids A. 4 (3): 633–636. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/1992PhFlA...4..633L">1992PhFlA...4..633L</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1063%2F1.858280">10.1063/1.858280</a>.</cite>

</li>

<li id="cite_note-Meneveauetal_1996-26"><a href="#cite_ref-Meneveauetal_1996_26-0">^</a> <cite id="CITEREFMeneveau,_C.Lund,_T._S.Cabot,_W._H.1996" class="citation journal cs1">Meneveau, C.; Lund, T. S.; Cabot, W. H. (1996). "A Lagrangian dynamic subgrid-scale model of turbulence". J. Fluid Mech. 319 (1): 353–385. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/1996JFM...319..353M">1996JFM...319..353M</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1017%2FS0022112096007379">10.1017/S0022112096007379</a>. <a href="Hdl_(identifier)" class="mw-redirect" title="Hdl (identifier)">hdl</a>:<a rel="nofollow" class="external text" href="https://hdl.handle.net/2060%2F19950014634">2060/19950014634</a>. <a href="S2CID_(identifier)" class="mw-redirect" title="S2CID (identifier)">S2CID</a> <a rel="nofollow" class="external text" href="https://api.semanticscholar.org/CorpusID:122183534">122183534</a>.</cite>

</li>

<li id="cite_note-Ghosal_1995-27"><a href="#cite_ref-Ghosal_1995_27-0">^</a> <cite id="CITEREFGhosalLundMoinAkselvoll1995" class="citation journal cs1">Ghosal, S.; Lund, T.S.; Moin, P.; Akselvoll, K. (1995). "A dynamic localization model for large-eddy simulation of turbulent flows". Journal of Fluid Mechanics. 286: 229–255. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/1995JFM...286..229G">1995JFM...286..229G</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1017%2FS0022112095000711">10.1017/S0022112095000711</a>. <a href="S2CID_(identifier)" class="mw-redirect" title="S2CID (identifier)">S2CID</a> <a rel="nofollow" class="external text" href="https://api.semanticscholar.org/CorpusID:124586994">124586994</a>.</cite>

</li>

<li id="cite_note-Carati_1995-28"><a href="#cite_ref-Carati_1995_28-0">^</a> <cite id="CITEREFCaratiGhosalMoin1995" class="citation journal cs1">Carati, D.; Ghosal, S.; Moin, P. (1995). <a rel="nofollow" class="external text" href="https://doi.org/10.1063%2F1.868585">"On the representation of backscatter in Dynamic Localization models"</a>. Physics of Fluids. 7 (3): 606–616. <a href="Bibcode_(identifier)" class="mw-redirect" title="Bibcode (identifier)">Bibcode</a>:<a rel="nofollow" class="external text" href="https://ui.adsabs.harvard.edu/abs/1995PhFl....7..606C">1995PhFl....7..606C</a>. <a href="Doi_(identifier)" class="mw-redirect" title="Doi (identifier)">doi</a>:<a rel="nofollow" class="external text" href="https://doi.org/10.1063%2F1.868585">10.1063/1.868585</a>.</cite>

</li>

</ol></div></div></div><div><div class="zim-footer">

This article is issued from <a class="external text" title="Last edited on 2025-03-05" href="https://en.wikipedia.org/wiki/?title=Large_eddy_simulation&oldid=1278965485">Wikipedia</a>. The text is available under <a class="external text" href="https://creativecommons.org/licenses/by-sa/4.0/deed.en">Creative Commons Attribution-Share Alike 4.0</a> unless otherwise noted. Additional terms may apply for the media files.

</div>

</div></div>

</div>

</main>

</div>

</body></html>